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ABSTRACT 

Entities  in  some  simulations  of  military  operations  move  unrealistically  from 
point  to  point  and  are  not  constrained  by  their  turning  radius.  The  fidelity 
of  this  representation  may  be  insufficient  for  operations  research  studies.  In 
this  paper  a  pursuer  intercepting  a  target  is  considered,  where  the  pursuer 
and  target  are  moving  at  constant  speeds  in  two  dimensions  and  the  target 
has  a  constant  velocity.  The  minimum  feasible  path  to  interception  for  a  given 
turning  radius  is  sought.  A  rigorous  analysis  of  the  model  constraints  produced 
an  algorithm  that  can  be  used  to  systematically  search  the  feasible  region  for 
the  minimum  path  to  interception.  At  the  core  of  the  algorithm  is  a  single 
implicit  equation  for  the  minimum  time  to  interception.  This  enables  the  effect 
of  turning  radius  to  be  incorporated  as  a  constraint  into  simulations  of  military 
operations,  improving  their  fidelity.  The  algorithm  is  also  straightforward  to 
implement  when  compared  with,  for  example,  a  traditional  flight  dynamics 
model,  and  has  a  broad  range  of  applications  in  path  optimisation  problems, 
the  development  of  computer  games  and  robotics. 
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Minimum  Paths  to  Interception  of  a  Moving  Target  when 
Constrained  by  Turning  Radius 

Executive  Summary 

Entities  in  some  simulations  of  military  operations  move  unrealistically  from  point  to  point 
and  are  not  constrained  by  their  turning  radius.  The  fidelity  of  this  representation  may 
be  insufficient  for  operations  research  studies. 

In  this  paper  a  pursuer  intercepting  a  target  is  considered,  where  the  pursuer  and 
target  are  moving  at  constant  speeds  in  two  dimensions  and  the  target  has  a  constant 
velocity.  The  minimum  feasible  path  to  interception  for  a  given  turning  radius  is  sought. 
The  pursuer  must  obey  three  constraints: 

(Cl)  The  pursuer  cannot  reverse  its  direction  to  intercept  a  target. 

(C2)  The  pursuer  cannot  intercept  a  target  inside  its  turning-circle. 

(C3)  The  pursuer  may  perform  at  most  one  complete  turn. 

It  is  assumed  that  the  pursuer’s  speed  is  strictly  greater  than  the  target’s  speed.  This 
assumption  is  not  absolutely  necessary,  however,  it  substantially  simplifies  the  analysis 
and  discussion. 

In  the  present  work,  a  rigorous  analysis  of  Constraints  (C1)-(C3)  produced  an  algo¬ 
rithm  that  can  be  used  to  systematically  search  the  feasible  region  for  the  minimum  path 
to  interception.  At  the  core  of  the  algorithm  is  a  single  implicit  equation  for  the  minimum 
time  to  interception.  This  equation  is  valid  in  an  arbitrary  Cartesian  coordinate  system 
and  encompasses  both  left  and  right  turns.  For  the  purpose  of  validation,  the  algorithm 
has  been  implemented  as  a  Mathematica  package  that  displays  the  minimum  feasible  path 
to  interception. 

Three  methods  are  proposed  for  incorporating  classification  range  into  the  present 
model:  an  exact  method,  a  heuristic  method,  and  a  method  that  includes  an  angle  of 
approach.  These  methods  provide  simple  models  of  a  pursuer’s  sensor  performance.  The 
present  model  can  be  easily  modified  to  encompass  the  heuristic  and  angle  of  approach 
methods. 

The  point-to-point  and  unconstrained  movement  of  entities  in  some  simulations  of 
military  operations  is  unrealistic.  The  algorithm  developed  here  enables  the  effect  of 
turning  radius  to  be  incorporated  as  a  constraint  into  these  simulations,  improving  their 
fidelity.  The  algorithm  is  also  straightforward  to  implement  when  compared  with,  for 
example,  a  traditional  flight  dynamics  model,  and  has  a  broad  range  of  applications  in 
path  optimisation  problems,  the  development  of  computer  games  and  robotics. 
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Notation 

i,  j,  k  standard  basis  vectors  in  the  x,  y  and  2:  directions 
X  a  generic  3-d  vector,  x  =  (xi,  X2,  X3) 

|x|  magnitude  of  x,  |x|  =  y^xf  -|-  x^  -|-  Xg 

X  •  y  scalar  (dot)  product,  x  •  y  =  xiyi  -|-  X2y2  +  x^ys 

X  X  y  vector  (cross)  product,  x  x  y  =  {x2y3  -  3:3^2,  -  Xiy3,xiy2  -  X2yi) 

n  n  =  0  results  in  a  left  turn,  n  =  1  results  in  a  right  turn 

Xc  centre  of  the  pursuer’s  turning-circle 

Tc  radius  of  the  pursuer’s  turning-circle 

Vci  pursuer’s  (dimensionless)  classification  range 

Xin  initial  position  of  the  pursuer 

Up  initial  velocity  of  the  pursuer 

Cp  speed  of  the  pursuer 

Xout  position  where  the  pursuer  exits  its  turn 

x!^  initial  position  of  the  target 

xt(0  position  of  the  target  at  time  t 

ut  velocity  of  the  target 

Ct  speed  of  the  target 

X/  point  of  interception 

Xci  point  of  classification 

a  |xl],  —  Xc|^  —  1  (in  dimensionless  form) 

P  ux  •  (xl!(,  —  Xc)  (in  dimensionless  form) 

e  Cx  /  C*p 

tc  time  for  the  pursuer  to  complete  its  turn 

ti  time  for  the  pursuer  to  move  from  Xout  to  x/ 

T  total  time  to  interception,  tc  +  M 

Tci  total  time  to  classification 

Tx  time  for  the  target  to  enter  the  pursuer’s  turning-circle  [see  Equation  (21)] 

Tfj  time  for  the  target  to  exit  the  pursuer’s  turning-circle  [see  Equation  (21)] 

To  time  to  interception  if  the  pursuer  does  not  turn  [see  Equation  (22)] 

To  time  to  interception  if  the  pursuer  does  not  turn  [see  Equation  (A5)] 

Ttt  time  to  interception  if  the  pursuer  takes  tt  units  of  time  to  turn  [see  Equation  (23)] 

Ttt  time  to  interception  if  the  pursuer  takes  tt  units  of  time  to  turn  [see  Equation  (24)] 

6c  angle  between  Xin  —  Xc  and  Xout  — 
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1  Introduction 


Entities  in  some  simulations  of  military  operations  move  unrealistically  from  point  to  point 
and  are  not  constrained  by  their  turning  radius.  The  fidelity  of  this  representation  may 
be  insufficient  for  operations  research  studies. 

In  this  paper  a  pursuer  intercepting  a  target  is  considered,  where  the  pursuer  and 
target  are  moving  at  constant  speeds  in  two  dimensions  and  the  target  has  a  constant 
velocity.  The  minimum  feasible  path  to  interception  for  a  given  turning  radius  is  sought. 
The  aim  is  to  rigorously  develop  an  algorithm  that  enables  the  effect  of  turning  radius  to 
be  incorporated  as  a  constraint  into  simulations  of  military  operations,  improving  their 
fidelity.  This  algorithm  should  be  straightforward  to  implement  when  compared  with,  for 
example,  a  traditional  flight  dynamics  model,  and  have  a  broad  range  of  applications  in 
path  optimisation  problems,  the  development  of  computer  games  and  robotics. 

During  maritime  surveillance  operations,  aircraft  search  areas  of  interest  in  order  to 
classify  as  many  ships  as  possible  in  the  shortest  possible  time.  Marlow,  Kilby  &  Mercer 
[2007]  endeavour  to  optimise  maritime  surveillance  operations  by  comparing  the  effect 
of  various  search  algorithms  on  the  surveillance  aircraft’s  performance.  At  present  their 
model’s  aircraft  simply  flies  from  point  to  point  and  is  not  constrained  by  its  turning 
radius. 

The  impact  of  turning  radius  and  classification  range  on  the  optimal  route  length  in 
a  simplified  maritime  surveillance  scenario  has  been  modelled  by  Mercer  et  al.  [2008], 
using  a  method  based  on  trigonometry.  A  similar  method  is  also  used  by  computer  game 
developers  [Pinter  2001,  Pinter  2002].  Both  of  these  methods  require  solving  a  system  of 
nonlinear  equations,  where  the  feasible  region  is  not  determined  and  the  resulting  curves 
are  not  necessarily  minimum  feasible  paths  to  interception. 

Classical  pursuit  curves  are  obtained  when  a  point  A  moves  with  constant  speed  to¬ 
wards  another  point  B  moving  with  constant  speed  along  a  known  curve,  where  the  tangent 
vector  at  A  is  required  to  be  parallel  to  the  line  connecting  A  and  B  [Boole  1859,  Colman 
1991,  Eliezer  &  Barton  1992,  Eliezer  &  Barton  1995,  Barton  &  Eliezer  2000,  Weisstein 
2008].  In  this  case,  the  turning  radius  of  A  is  not  constrained  and  the  resulting  curves 
are  also  not  necessarily  minimum  feasible  paths  to  interception.  In  the  present  work,  the 
tangent  vector  of  the  pursuer  is  not  required  to  be  parallel  to  the  line  connecting  the 
pursuer  to  the  target. 

The  objective  of  pursuit-evasion  games  is  to  determine  the  optimal  strategies  that 
result  in  a  pursuer  capturing  an  evader.  The  theory  of  pursuit-evasion  games  was  first 
studied  by  Isaacs  [1965]  and  has  been  applied  to  a  variety  of  problems  in  the  guidance  and 
optimal  control  literature:  for  example,  see  Shinar,  Guelman  &  Green  [1989]  and  Shima 
&  Shinar  [2002].  The  present  work  was  motivated  by  research  into  maritime  surveillance 
operations  [Marlow,  Kilby  &  Mercer  2007] .  In  this  instance,  the  speed  of  the  surveillance 
aircraft  is  typically  much  greater  than  the  speed  of  the  ship,  and  consequently  any  evasive 
manoeuvres  executed  by  the  ship  will  have  little  impact  on  the  surveillance  aircraft’s  ability 
to  classify  the  ship.  Eor  this  reason,  the  target  considered  here  does  not  attempt  to  evade 
the  pursuer,  and  although  this  can  be  regarded  as  a  special  case  of  a  pursuit-evasion  game, 
a  simpler,  more  direct  method  will  be  used  to  construct  a  solution. 
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Figure  1:  Definition  diagram  showing  two  possible  paths  to  intereeption  resulting  from 
left  and  right  turns.  The  pursuer  (P)  turns  to  intereept  the  target  (T)  at  time  t  =  t).  The 
pursuer’s  initial  veloeity  is  Up,  and  the  target’s  veloeity  is  ux-  The  pursuer  and  target 
have  eonstant  speeds  Cp  and  Ct,  respeetively. 


Dubins  [1957]  studied  planar  continuously  differentiable  curves  of  minimal  length  with 
average  curvature  bounded  by  R~^,  between  prescribed  initial  and  final  positions  and 
orientations.  He  proved  that  such  curves  exist  and  are  necessarily  a  subpath  of  a  path 
of  type  CLC  or  of  type  CCC,  where  C  is  an  arc  of  a  circle  of  radius  i?,  and  L  is  a 
straight  line  segment.  These  results  can  be  applied  to  the  motion  planning  of  a  car-like 
robot  that  only  moves  forwards.  Reeds  &:  Shepp  [1990]  have  extended  the  work  of  Dubins 
[1957]  to  allow  for  both  forwards  and  backwards  motion  (that  is,  paths  with  cusps),  and 
provide  explicit  formulae  for  the  resulting  68  paths.  There  are  many  other  generalisations 
of  Dubins’  paper:  for  example,  to  the  motion  planning  of  a  car-like  robot  moving  amid 
obstacles  [Laumond  et  al.  1994],  and  to  the  Traveling  Salesperson  Problem  for  a  car¬ 
like  robot  [Savla,  Frazzoli  &:  Bullo  2008].  The  paths  considered  in  this  paper  differ  from 
Dubins’  paths,  as  the  pursuer’s  final  position  and  orientation  are  not  prescribed. 
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Not  all  paths  to  interception  will  be  considered  in  this  paper.  Instead,  as  the  final 
orientation  of  the  pursuer  is  not  prescribed  and  since  the  shortest  path  between  two  points 
is  a  straight  line,  the  work  of  Dubins  [1957]  suggests  that  minimal  paths  to  interception 
will  consist  of  an  arc  of  a  circle  of  radius  equal  to  the  turning  radius  of  the  pursuer,  followed 
by  a  straight  line  segment. 

The  definition  diagram  for  the  present  scenario  is  shown  in  Figure  1.  At  time  t  =  0  the 
pursuer  at  position  Xin  with  velocity  Up  turns  with  a  turning  radius  of  Tc  to  intercept  the 
target,  which  is  initially  at  position  xlj.  with  a  constant  velocity  ux.  It  may  be  possible 
for  the  pursuer  to  turn  either  left  or  right  to  intercept  the  target.  The  pursuer  exits 
its  turn  at  position  Xout  {t  =  tc)  and  intercepts  the  target  at  position  x/  {t  =  tc  +  ti). 
The  pursuer  and  target  move  in  two  dimensions  and  have  constant  speeds  Cp  and  Cx, 
respectively.  The  minimum  feasible  time  to  interception,  and  hence  the  minimum  feasible 
path  to  interception,  for  a  given  turning  radius  is  sought.  The  pursuer  must  obey  three 
constraints: 

(Cl)  The  pursuer  cannot  reverse  its  direction  to  intercept  a  target. 

(C2)  The  pursuer  cannot  intercept  a  target  inside  its  turning-circle. 

(C3)  The  pursuer  may  perform  at  most  one  complete  turn. 

It  is  assumed  that  the  pursuer’s  speed  is  strictly  greater  than  the  target’s  speed.  This 
assumption  is  not  absolutely  necessary,  however,  it  substantially  simplifies  the  analysis 
and  discussion. 

In  maritime  surveillance  operations,  aircrew  can  classify  a  ship  if  the  aircraft  is  within 
a  certain  distance  of  the  ship,  which  depends  on  the  aircraft’s  sensor  suite,  weather,  sea 
state,  et  cetera.  Three  methods  will  be  proposed  for  incorporating  classification  range 
into  the  present  scenario:  an  exact  method,  a  heuristic  method,  and  a  method  that  in¬ 
cludes  an  angle  of  approach.  These  methods  provide  simple  models  of  a  pursuer’s  sensor 
performance. 
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2  The  mathematical  model 

In  this  section,  a  mathematical  model  of  the  scenario  described  in  the  Introduction  is 
derived. 

2.1  Derivation  of  the  governing  equations 

Any  point  x  on  a  circle  of  radius  Tc  centred  on  Xc  satisfies 

|x  -  Xc|  =  re- 

Since  Xc  —  Xin  is  orthogonal  to  Up  (see  Figure  1)  it  follows  that  Xc  satisfies 

|xin  -  Xc|  =  Tc, 

(Xc  -  Xin)  •  Up  =  0. 

Using  the  aforementioned  reasoning,  the  exit  point  on  the  turning-circle  obeys 

l^out  Xc|  — 

(^C  ^out)  ■  (x/  ^out)  -  0, 

where  the  point  of  interception  is  given  by 

X/  =  xJ]^  -|-  T ux  • 

Here  T  =  tc  +  ti  is  the  total  time  to  interception. 

Since  the  pursuer  is  moving  at  a  constant  speed,  the  time  taken  for  the  pursuer  to 
complete  the  turn  is 

_  OcTc 
“  Cp  ’ 

where  9c  is  the  angle  between  Xin  —  Xc  and  Xout  ~  (see  Figure  1).  The  definition  of  the 
scalar  product  yields  an  expression  for  6c  [Spiegel  1974]: 

\  (xin  Xc)  •  (Xout  Xc) 

cos(6'c)  = - 5 - . 

ri 

Once  again  since  the  pursuer  is  moving  at  a  constant  speed,  the  time  taken  for  the  pursuer 
to  move  from  the  exit  point  to  the  point  of  interception  is 

|x/-Xout| 

For  convenience  and  to  highlight  the  key  parameters  of  the  problem,  we  introduce  the 
following  dimensionless  variables: 

X  =  TeX,  t  =  ^i,  Up  =  CpUp,  Ux  =  CxUx,  (1) 

Op 
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where  a  caret  indicates  a  dimensionless  variable  and  x  is  a  generic  position  vector.  Em¬ 
ploying  these  scales  results  in  the  following  dimensionless  system: 


Xin  -  Xc  =  1, 

(2) 

(Xc  -  Xin)  •  =  0, 

(3) 

l^out  ^c]  -  1) 

(4) 

(Xc  -  Xout)  •  (x/  -  Xout)  =  0, 

(5) 

X/  =  x^  -b  cTut, 

(6) 

T  =  tc  +  ti, 

(7) 

cos(fc)  =  (xin  -  Xc)  •  (Xout  “  Xc), 

(8) 

tl  —  jx/  Xout]  ) 

(9) 

where  the  dimensionless  parameter  e  is  defined  by^ 

(7x 

Note  that  the  carets  have  been  omitted  for  convenience,  thus  all  functions,  variables  and 
parameters  will  henceforth  refer  to  dimensionless  quantities. 


2.2  Solution  of  the  governing  equations 

Vectors  of  the  form  X  =  (a,  b,  0)  have  the  following  properties: 

1.  |X  X  k|  =  |X|; 

2.  (X  X  k)  -  X  =  0; 

3.  (X  X  k)  X  k  =  -X; 

4.  (X  X  k)  X  X  =  |X|^k, 

where  k  =  (0,0, 1).  Properties  1  and  2  can  be  verified  by  inspection.  Properties  3  and  4 
can  be  derived  using  the  following  identity  for  general  3-dimensional  vectors  [Spiegel  1974]: 

(A  X  B)  X  C  =  (A  •  C)  B  -  (B  •  C)  A. 

Properties  1  to  4  will  now  be  used  to  solve  Equations  (2)  to  (9). 

stationary  target  implies  that  e  =  0.  In  this  case  Equations  (2)  to  (9)  decouple  and  an  explicit 
solution  for  the  total  time  to  interception  can  be  obtained.  An  approximate  expression  has  been  derived 
that  is  valid  in  the  limit  e  ^  0.  However,  this  expression  is  cumbersome  and  only  accurate  for  very  small 
values  of  e. 


5 


DSTO-TR-2227 


2.2.1  The  exit  point  on  the  tnrning-circle 


It  can  be  verified  using  Properties  1  and  2  that  the  solution  to  Equations  (2)  and  (3)  is 

Xc  =  Xin  +  (-l)”+^u^  X  k,  (10) 

as  I  Up  I  =  1,  where  n  =  0  or  n  =  1  and  determines  the  direction  of  the  turn.^  Similarly, 
an  implicit  solution  to  Equations  (4)  and  (5)  is 


(-1)” 

Xout  —  Xc  +  -  (x/  Xout)  X  k, 


(11) 


where  Equation  (9)  has  been  used.  Next,  Property  3  and  Equation  (11)  give 


Xout  X  k  =  Xc  X  k  - 


(-1)" 

ti 


(X/  -  Xout)  • 


Substituting  this  result  into  Equation  (11)  and  solving  for  Xout  leads  to  an  explicit  general 
solution  for  the  exit  point: 


^OUt 


^  2  +  (“1)”^;  -  ^c)  X  k) . 

1  +  t; 


(12) 


2.2.2  The  total  time  to  interception 


Equation  (9)  and  Pythagoras’  Theorem  yield  (see  Eigure  1) 


tl  =  Y^|x/  -  Xc|^  -  1. 

(13) 

This  result  in  conjunction  with  Equation  (6)  leads  to 

ti{T)  =  V(eT)2  +  2PeT  +  a. 

(14) 

as  ut  =  1,  where 

a  =  |xt  —  Xc|^  —  1, 

(15) 

/3  =  Ut  •  (x)].  -  Xc)  . 

(16) 

Note  that  the  target  is  initially  outside  the  pursuer’s  turning-circle  if  and  only  if  a  >  0. 
The  parameter  —j3  is  the  component  of  the  target’s  heading  in  the  direction  of  Xc,  which 
follows  from  the  definition  of  the  scalar  product. 

Equations  (12)  and  (14)  reveal  that 

Xout(r)  =  (x/(T)  +  tz2(r)xe  +  MT)  -  X,)  X  k)  .  (17) 

Observe  that  the  only  unknown  in  Xout  is  T.  As  a  result  the  expression  for  tc  [Equation  (8)] 
can  be  recast  as 

COs(r  -  tl{T))  =  (Xin  -  Xc)  •  (Xout(T)  -  Xc),  (18) 

since  tc  =  T  —  ti.  Hence  the  governing  system  [Equations  (2)  to  (9)]  has  been  reduced 
to  a  single  implicit  equation  for  the  total  time  to  interception.  Equation  (18)  is  valid  in 
an  arbitrary  Cartesian  coordinate  system  and  encompasses  both  left  and  right  turns.  The 
minimal  solution  to  Equation  (18)  is  sought  such  that  Constraints  (C1)-(C3)  are  satisfied. 

n  =  0  then  the  pursuer  will  turn  left,  whereas  if  n  =  1  then  the  pursuer  will  turn  right. 
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2.3  Derivation  of  the  constraints 


In  this  section,  mathematical  representations  of  Constraints  (C1)-(C3)  are  derived,  and 
the  impact  of  these  constraints  on  the  feasibility  of  paths  to  interception  is  considered 
one-at-a-time. 

(Cl):  The  pursuer  cannot  reverse  its  direction  on  the  turning-circle  to  intercept  a 
target,  which  can  be  expressed  mathematically  as 


sgn((xc  -  Xin)  X  •  k)  =  sgn((xc  -  Xout)  x  (x/  -  Xout)  •  k)  ,  (19) 


where 


sgn(a;) 


—  1,  X  <  0 

^  0,  X  =  0 

1,  X  >  0. 

\  ’ 


Equation  (19)  states  that  the  centre  of  the  turning-circle  must  remain  on  the  same  side  of 
the  pursuer  as  it  enters  and  exits  the  turning-circle. 


Using  Property  4  from  Section  2.2  together  with  Equations  (11)  and  (13),  it  can  be 
shown  that 


(Xc  -  Xin)  X  •  k  =  (-1)''+^ 

(^c  ^out)  X  (x/  Xout)  *  k  —  (  1)  tl^ 

and  hence  Xout  satisfies  Constraint  (Cl)  by  construction.  Furthermore,  the  turning-circle 
centred  on  Xc  with  n  =  0  has  an  exit  point  with  n  =  0,  and  the  other  turning-circle  centred 
on  Xc  with  n  =  1  has  an  exit  point  with  n  =  1. 

(C2):  The  pursuer  cannot  intercept  the  target  inside  its  turning-circle,  that  is, 

|X/  -  Xc|  ^  1, 


which  evaluates  to 

{eTf +  2peT  +  0,  (20) 

that  is,  ti{T)  must  be  a  real  function  [see  Equation  (14)].  The  zeros  of  ti{T)  are 


Tl 


e 


Tr 


-(3  -I-  (3“^  -  a 

e 


(21) 


The  zeros  Tl  and  Tr  are  both  finite,  real  and  positive  when  a  >  0  and  e  >  0  and  0^  a 
and  /3  <  0,  which  follows  from  \(3\  >  0 0^  —  a. 

It  can  be  seen  from  Equation  (20)  that  Constraint  (C2)  is  satisfied  when  a  >  0  and 
e  =  0,  which  corresponds  to  a  stationary  target  that  is  initially  outside  the  pursuer’s 
turning-circle.  Now  let  a  >  0  and  0  <  e  <  1.  In  this  case  Equations  (20)  and  (21)  reveal 
Constraint  (C2)  is  satisfied  when  /3  ^  0  or  ^  a.  The  condition  ^  0  implies  the 
target  is  heading  away  from  the  centre  of  the  turning  circle,  which  follows  directly  from 
the  definition  of  (3.  If  ^  a  is  true  then  the  target  will  never  enter  the  interior  of  the 
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pursuer’s  turning-circle,  which  can  be  verified  as  follows.  Let  XT(t)  =  xl],  -|-  etuT  be  the 
position  of  the  target  at  time  t.  The  target  will  never  enter  the  interior  of  the  pursuer’s 
turning-circle  if  and  only  if  |xT(t)  —  Xc|  ^  1  for  all  t.  In  this  instance  Equation  (20)  will 
be  satisfied  and  hence  (eT)^  -|-  2f5eT  +  a  will  have  at  most  one  real  zero,  that  is,  ^  a. 
Therefore,  if  a  >  0  and  0  <  e  <  1  and  /3  <  0  and  >  a,  the  target  will  enter  the  pursuer’s 
turning-circle  at  time  Tl  and  exit  the  turning-circle  at  time  Tr. 

To  summarise.  Constraint  (C2)  is  satisfied  if  the  target  is  initially  outside  the  pursuer’s 
turning-circle  and  will  never  enter  the  interior  of  the  turning-circle  (a  >  0  and  e  =  0  or 
^  0  or  ^  a).  If  the  target  is  initially  outside  the  pursuer’s  turning-circle  and  will 
enter  the  interior  of  the  turning-circle  at  some  time  (a  >  0  and  0  <  e  <  1  and  /3  <  0  and 

>  a),  then  a  solution  to  Equation  (18)  that  satisfies  Constraint  (C2)  may  exist  in  the 
domain  0  ^  T  ^  Ti  or  T  ^  Tr. 

A  discussion  of  the  case  where  the  target  is  initially  inside  (or  on  the  boundary  of)  the 
pursuer’s  turning-circle  (a  ^  0)  is  postponed  until  Appendix  A. 3. 

(C3):  The  pursuer  may  perform  at  most  one  complete  turn,  implying  0  ^  tc  ^  27r. 
However  it  is  only  necessary  to  search  for  solutions  in  0  ^  tc  ^  tt,  which  can  be  shown 
as  follows.  Given  that  Xout  is  defined  on  a  circle,  it  is  a  27r-periodic  function  of  tc,  as  is 
cosine.  Therefore  if  vr  ^  tc  ^  27r  is  a  solution  of  Equation  (8)  then  so  is  0  ^  27r  —  tc  ^  tt, 
implying  the  minimal  solution  will  occur  in  the  interval  0  ^  tc  ^  vr.^ 

Given  the  aforementioned  discussion  and  since  tc(T)  =  T  —  ti{T),  a  solution  to  Equa¬ 
tion  (18)  must  be  restricted  to  the  interval  Tq  ^  T  ^  T^r,  where 


e/3-h  V(e/3)2 -h  (1  -  e2)a 

(22) 

e(5  +  Ti  +  ^ -  a) 

-  1  _  e2 

(23) 

_  efd +  Tt  -  +  €^{(3“^  -  a) 

^TT  1  9  ’ 

1  — 

(24) 

which  are  obtained  by  solving 

1 

§ 

II 

(25) 

1 

II 

(26) 

-  ti(T^)  =  TT. 

(27) 

Therefore  Tq  is  the  time  to  interception  that  corresponds  to  the  pursuer  not  turning. 
Likewise,  and  are  the  times  to  interception  that  correspond  to  the  pursuer  taking 
TT  units  of  time  to  turn.  Refer  to  Appendix  A  for  a  discussion  of  the  validity  of  these 
solutions  of  Equations  (25)  to  (27). 

®An  exception  occurs  if  0  ^  tc{T)  rr  is  the  minimal  solution  of  Equation  (8)  and 

(xout(T)  —  Xin)  •  Up  ^  0,  then  the  minimum  time  to  complete  the  turn  is  27r  —  tc-  This  follows  from 
Figure  1  and  the  dehnition  of  the  scalar  product;  see  Section  4.1.2  for  more  details. 
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If  To  and  are  real  numbers,  then  Tq  ^  T  ^  T^-  is  a  necessary  condition  to  satisfy 
Constraint  (C3),  however  it  may  not  be  sufficient.  Since  tc{T)  =  T  —  ti{T),  if  Con¬ 
straint  (C2)  is  not  satisfied  then  Constraint  (C3)  cannot  be  satisfied.  Further  discussion 
of  the  feasibility  of  times  to  interception  with  respect  to  Constraint  (C3)  is  postponed 
until  Section  3. 

It  can  be  seen  from  Equations  (22)  to  (24)  that  there  is  a  singularity  at  e  =  1.  Although 
solutions  of  Equations  (2)  to  (9)  exist  for  the  unusual  case  of  e  ^  1,  these  solutions  have 
a  different  form  to  those  obtained  for  0  ^  e  <  1.  Since  the  case  0  ^  e  <  1  is  of  greater 
practical  interest,  for  the  reminder  of  the  paper  it  is  assumed  that  0  ^  e  <  1. 


9 


DSTO-TR-2227 


3  Feasible  times  to  interception 

Feasible  times  to  interception  are  defined  to  be  those  times  which  simultaneously  satisfy 
Constraints  (C1)-(C3).  However,  there  is  room  for  confusion  here  as  Equation  (18)  is 
both  an  additional  constraint  and  the  objective  function.  This  is  because  a  solution  to 
Equation  (18)  may  not  exist  from  the  set  of  feasible  times  to  interception  (called  the 
feasible  region),  and  hence  Equation  (18)  acts  as  an  additional  constraint.  If  solutions  to 
Equation  (18)  do  exist  from  the  feasible  region,  then  Equation  (18)  becomes  the  objective 
function  and  the  minimum  solution  is  chosen. 

A  rigorous  analysis  of  Constraints  (C1)“(C3)  has  been  performed  in  Appendix  A.  The 
results  from  this  appendix  and  Section  2.3  will  now  be  utilised  to  determine  the  feasible 
region.  To  simplify  the  discussion,  only  the  case  where  the  target  is  initially  outside  the 
pursuer’s  turning-circle,  that  is,  a  >  0,  will  be  considered  here;  the  case  when  a  ^  0  is 
analyzed  in  Appendix  A. 3. 

In  Section  2.3  it  was  shown  that  Constraint  (Cl)  is  satisfied  by  construction,  and  (if 
a  >  0)  there  are  two  main  cases  effecting  the  feasibility  of  times  to  interception  with 
respect  to  Constraints  (C2)  and  (C3).  These  cases  are: 

•  the  target  will  never  enter  the  interior  of  the  pursuer’s  turning-circle  (a  >  0  and 
e  =  0  or  /3  ^  0  or  ^  a),  and; 

•  the  target  will  enter  the  interior  of  the  pursuer’s  turning-circle  at  some  time  (a  >  0 
and  0  <  e  <  1  and  /3  <  0  and  >  a). 

Hence  the  parameters  effecting  the  feasibility  of  times  to  interception  are  a  =  |  —  Xc  |  ^  —  1 , 

/?  =  ux  •  (xlj.  —  Xc)  and  e  =  Cx/Cp.  Eurthermore,  when  a  >  0,  it  is  natural  to  partition 
the  parameter  space  into  the  disjoint  sets 

{ (a,  /3,  e)  I  e  =  0  or  /3  ^  0  or  ^  a}  , 


and 

{(a,  /?,  e)  I  0  <  e  <  1  and  /3  <  0  and  0^  >  a]  . 

Proposition  7  of  Appendix  A. 2  states  that  if  a  >  0  and  e  =  0  or  /3  ^  0  or  ^  a, 
then  feasible  times  to  interception  T  occur  in  the  interval  Tq  ^  T  ^  T^r.  In  this  instance, 
the  target  will  never  enter  the  interior  of  the  pursuer’s  turning-circle  and  so  the  pursuer 
is  free  to  perform  a  complete  turn  before  attempting  an  intercept.  In  fact,  it  is  proven 
in  Proposition  8  of  Appendix  A. 2  that  the  pursuer  will  have  at  least  one  opportunity  to 
intercept  the  target. 

The  other  case,  when  a  >  0  and  0  <  e  <  1  and  /?  <  0  and  >  a,  is  more  complex.  In 
this  situation  there  are  four  conditions  effecting  the  feasibility  of  times  to  interception  that 
lead  to  four  different  feasible  regions,  which  are  presented  in  Table  1.  The  results  in  this 
table  are  proven  in  Propositions  9  to  12  of  Appendix  A. 2.  These  propositions  establish 
the  shape  of  tc{T)  which  is  shown  in  Eigure  2(a)-(d),  where  the  four  feasible  regions  in 
Table  1  can  also  be  obtained  by  inspection.  Recall  that  in  all  four  cases  the  target  will 
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Case 

Condition 

Eeasible  Region 

1 

Tr^tt 

To  ^T  ^Tr  and  Trs^T^T^, 

2 

Tl  <  TT  <  Tr  and  is  real 

To  ^Tr  and  ^  T  ^ 

3 

Tr  <  TT  <  Tr  and  T^r  is  complex 

To  ^T  if  Tr 

4 

TlP  TT 

To  ^T  ^Tt^ 

Table  1:  The  four  cases  effecting  the  feasibility  of  times  to  interception  when  the  target 
will  enter  the  interior  of  the  pursuer’s  turning- circle  at  some  time  (a  >  and  0  <  e  <  1 
and  P  <  0  and  >  a). 


enter  the  interior  of  the  pursuer’s  turning-circle  at  some  time.  Physical  meaning  can  then 
be  attributed  to  the  four  cases  as  follows:^ 

Case  1:  The  target  will  exit  the  pursuer’s  turning-circle  before  time  t  =  27r,  and  the 
pursuer  can  feasibly  intercept  the  target  after  it  exits  the  turning-circle.  Equivalently, 
the  target  is  moving  fast  enough  for  the  pursuer  to  attempt  an  intercept  immediately 
after  the  target  leaves  the  turning-circle. 

Case  2:  The  target  will  exit  the  pursuer’s  turning-circle  before  t  =  27r,  however  the 
pursuer  cannot  feasibly  intercept  the  target  after  it  exits  the  turning-circle  until 
t  =  Tn-  Equivalently,  the  target  is  moving  fast  enough  for  the  pursuer  to  attempt 
an  intercept  T^^  —  Tr  units  of  time  after  the  target  leaves  the  turning-circle. 

Case  3:  The  target  will  not  exit  the  pursuer’s  turning-circle  before  t  =  27r,  and  hence 
the  pursuer  cannot  feasibly  intercept  the  target  after  it  exits  the  turning-circle. 
Equivalently,  the  target  has  enough  speed  to  enter  the  pursuer’s  turning-circle,  but 
is  moving  too  slowly  to  leave  the  turning-circle  before  the  pursuer  has  performed  a 
complete  turn. 

Case  4:  The  target  will  not  enter  the  pursuer’s  turning-circle  before  t  =  27r,  and  hence 
the  pursuer  can  feasibly  perform  a  complete  turn  before  intercepting  the  target. 
Equivalently,  the  target  is  moving  so  slowly,  and/or  it  is  so  far  away  from  the  pur¬ 
suer’s  turning-circle,  that  it  will  not  enter  the  turning-circle  before  the  pursuer  has 
performed  a  complete  turn. 

The  above  cases  only  describe  the  feasibility  of  times  to  interception  with  respect  to 
Constraints  (C1)-(C3),  and  feasibility  does  not  imply  that  the  pursuer  can  intercept  the 
target.  However,  it  is  proven  in  Proposition  9  of  Appendix  A. 2  that  if  Tr  ^  tt  (Case  4), 
then  the  pursuer  will  have  at  least  one  opportunity  to  intercept  the  target. 


^The  speed  of  the  target  is  relative  to  the  speed  of  the  pursuer. 
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UT) 


(a)  If  Tr  ^  TT,  then  the  feasible  region  is  Tq  ^  T  ^  Tl  and  Tr  <  T  ^  T^.  (See  Proposition  10  in 
Appendix  A. 2  for  a  proof  of  this  statement.) 

tc{T) 


(b)  If  Tl  <  TT  <  Tr  and  is  real,  then  the  feasible  region  is  Tq  ^  T  <  Tr  and  ^  (See 

Proposition  12  in  Appendix  A. 2  for  a  proof  of  this  statement.) 

Figure  2:  The  time  for  the  pursuer  to  complete  its  turn  tc{T)  =  T  —  ti{T)  versus  the  total 
12  time  to  interception  T,  for  a  >  0  and  0  <  e  <  1  and  /?  <  0  and  >  a  (Cases  1  and  2). 
The  part  of  the  graph  that  is  displayed  as  a  solid  line  is  the  feasible  region. 
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UT) 


(c)  If  Tl  <  tt  <T[i  and  is  complex,  then  the  feasible  region  is  Tg  <  T  ^  T^.  (See  Proposition  11 
in  Appendix  A. 2  for  a  proof  of  this  statement.) 

tc{T) 


(d)  If  T]^  ^  TT,  then  the  feasible  region  is  Tq  ^  T  ^  (See  Proposition  9  in  Appendix  A. 2  for  a 
proof  of  this  statement.) 

Figure  2:  The  time  for  the  pursuer  to  eomplete  its  turn  tc(T)  =  T  —  ti(T)  versus  the  total 
time  to  intereeption  T,  for  a  >  0  and  0  <  e  <  1  and  /?  <  0  and  ff^  >  a  (Cases  3  and  f). 
The  part  of  the  graph  that  is  displayed  as  a  solid  line  is  the  feasible  region. 
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4  The  Interception  algorithm 

The  results  obtained  in  Section  3  have  been  utilised  to  develop  an  algorithm  that  can  be 
used  to  systematically  search  the  feasible  region  for  the  minimum  time  to  interception. 
The  algorithm,  called  Interception,  is  shown  in  Figure  3.  A  description  of  Interception 
using  words,  rather  than  symbols,  is  presented  in  Appendix  B. 

To  simplify  the  discussion,  only  the  case  where  the  target  is  initially  outside  the  pur¬ 
suer’s  turning-circle,  that  is,  a  >  0,  will  be  considered  in  this  section;  the  Interception 
algorithm  for  the  case  when  a  ^  0  is  derived  in  Appendix  A. 3  and  displayed  in  Fig¬ 
ure  Cl(a)-(b). 

4.1  Implementing  Interception 

In  this  section,  the  practicalities  of  implementing  Interception  are  discussed. 

4.1.1  Convergence  of  the  root-finding  method 

The  Interception  algorithm  can  be  used  to  systematically  search  the  feasible  region  for  the 
minimum  time  to  interception,  that  is,  the  minimum  feasible  solution  of  Equation  (18). 
However,  the  reliability  of  Interception  to  converge  to  the  minimum  time  to  interception 
is  limited  by  the  root-finding  method  used  to  solve  Equation  (18).  Eor  instance,  the 
root-finding  method  may  fail  to  converge,  even  when  a  solution  does  exist,  or  it  may  not 
converge  to  the  minimum  solution.  These  issues  can  be  alleviated  as  follows. 

Suppose  a  solution  to  Equation  (18)  is  sought  in  the  interval  L  ^  T  ^  R.  Then  choose 
(5i  and  f  such  that  di  >  1  and  0  <  (/>  <  1.  Since  the  minimum  solution  is  sought,  begin 
the  search  aX  L  +  {R  —  L) / 6i.  If  the  root-finding  method  fails  to  converge  to  a  solution 
in  the  interval  L  ^  T  ^  R,  then  set  62  =  (fSi  and  search  again,  this  time  starting  at 
L  +  {R  —  L) / 62-  Continue  in  this  manner,  starting  at  L  +  {R  —  L) / 5^  for  the  A^th  attempt 
{5n  =  until  the  root-finding  method  converges  to  a  solution  of  Equation  (18) 

in  the  interval  L  ^  T  ^  R  or  5m  ^  I-  If  a  solution  has  still  not  been  found,  depending 
on  which  root-finding  method  is  being  used,  it  may  be  necessary  to  repeat  this  procedure 
starting  the  search  at  R  —  {R  —  L) / 5i.  By  changing  the  parameters  di  and  cj)  the  balance 
between  reliability  and  computational  performance  may  be  adjusted. 

Root-finding  methods  for  nonlinear  equations  will  not  be  discussed  further  in  this 
paper  as  there  are  a  multitude  of  text  books  on  numerical  analysis  that  cover  the  topic; 
for  example,  see  Atkinson  [1989]. 

4.1.2  The  minimum  feasible  time  to  interception 

Although  Constraint  (C3)  implies  that  0  ^  tc  ^  27r,  it  is  only  necessary  to  search  for 
a  solution  to  Equation  (18)  in  0  ^  tc  ^  tt,  as  discussed  in  Section  2.3.  However,  the 
resulting  time  to  interception  will  not  equal  the  true  minimum  time  to  interception  if 
the  pursuer  performs  a  turn  of  more  than  tt  radians.  This  does  not  cause  any  problems. 
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Figure  3:  The  Interception  algorithm  for  determining  feasible  times  to  intereeption  T; 
the  ease  when  a  ^  0  zs  eontinued  on  Figure  Cl  (a).  A  deseription  of  Interception  using 
words,  rather  than  symbols,  is  presented  in  Appendix  B. 
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because  the  true  minimum  time  to  interception  can  be  easily  obtained  from  the  solution 
of  Equation  (18)  as  follows. 

The  true  minimum  time  for  the  pursuer  to  complete  its  turn  is 


T-ti{T),  (Xout(T)  -  Xin)  •  >  0 

27r  -  (T  -  ti{f)),  (xout(T)  -  Xin)  •  ^  0, 


(28) 


which  follows  from  Figure  1  and  the  definition  of  the  scalar  product.  The  true  minimum 
feasible  time  to  interception  is  then  given  by 


T  =  tc{T)+ti{T),  (29) 

where  T  is  the  minimum  feasible  solution  of  Equation  (18)  returned  by  Interception;  see 
Figure  3. 


4.1.3  Errors  and  warnings 


When  Interception  is  executed,  a  number  of  error  and  warning  messages  may  be  returned 
that  correspond  to  physical  events.^  These  messages  are  described  below: 


Error  1  is  returned  if 

•  e<0ore^l:  e<0is  not  physically  possible  and  e  ^  1  implies  the  target’s 
speed  is  greater  than  or  equal  to  the  speed  of  the  pursuer,  or; 

•  |up|  7^  1  or  |ux|  7^  1:  the  headings  have  been  scaled  such  that  |up|  =  |ux|  =  1.® 

Error  2  is  returned  if  the  root-finding  method  used  to  solve  Equation  (18)  fails  to  con¬ 
verge,  when  a  solution  is  known  to  exist;  see  Propositions  8  and  9  from  Appendix  A. 2. 

Warning  1  is  returned  if  the  root-finding  method  used  to  solve  Equation  (18)  fails  to 
converge.  In  this  instance,  the  pursuer  is  unable  to  intercept  the  target  before  it 
enters  the  turning-circle  (if  a  >  0),  and  the  target  is  moving  too  slowly  to  leave  the 
turning-circle  before  the  pursuer  has  performed  a  complete  turn;  see  Case  3  from 
Section  3. 

Warning  2  is  returned  if  the  root-finding  method  used  to  solve  Equation  (18)  fails  to 
converge.  In  this  instance,  the  pursuer  is  unable  to  intercept  the  target  before  it 
enters  the  turning-circle  (if  a  >  0),  and  the  target  is  moving  too  slowly  to  leave  the 
turning-circle  in  time  for  the  pursuer  to  attempt  an  intercept;  see  Cases  1  and  2 
from  Section  3. 


®Error  2  does  not  correspond  to  a  physical  event. 

®This  error  is  irrelevant  if  Interception  is  implemented  using  dimensional  quantities. 
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4.1.4  A  modified  definition  of  a 

The  overall  behaviour  of  Interception  is  determined  by  the  sign  of  a  (see  Figure  3).  It  has 
been  found  that  Interception  is  very  sensitive  to  numerical  error  if  |xS].  —  Xc|  ~  1,  which 
stems  from  the  definition  of  a: 

a  =  |x!^  —  Xc|^  —  1  =  (|xx  —  Xc|  —  l)  (|xx  —  Xc|  +  l)  . 

It  can  be  demonstrated  that  the  squared  term  introduces  an  additional  numerical  error  if 
|xl(i  —  Xc|  ~  1,  that  is,  if  |x!]^  —  Xc|  numerically  equals  unity,  a  will  not  numerically  equal 
zero.  The  following  modified  definition  of  a  alleviates  this  problem: 

/ 

I  Xc  I  1 ,  II  X.J.  Xc  I  f  I  ^  ^tol 

0 ,  II  Xrj.  Xc  I  f  I  ^  ^tol 

I  I  J-  I  I 

where  0  <  5toi  1  is  a  numerical  tolerance. 

4.1.5  Polar  coordinates 

It  may  be  useful  to  express  the  points  Xin  and  Xout  in  polar  coordinates  centred  on  Xc. 
A  difficulty  here  is,  standard  trigonometric  methods  for  transforming  points  in  Cartesian 
coordinates  into  polar  coordinates,  may  return  points  in  the  wrong  quadrants.  A  method 
to  obtain  expressions  for  Xin  and  Xout  in  polar  coordinates  centred  on  Xc,  such  that  the 
resulting  points  are  in  the  correct  quadrants,  is  presented  in  Appendix  E. 

4.2  An  implementation  of  Interception 

For  the  purpose  of  validation,  the  Interception  algorithm  has  been  implemented  as  a 
Mathematical  package,  called  TurningCircle  .m,  that  displays  the  minimum  feasible  path 
to  interception.  The  code  for  TurningCircle  .m  is  included  in  Appendix  D. 

Whether  or  not  TurningCircle  .m  returns  the  minimum  path  to  interception  depends 
on  the  reliability  of  the  root-finding  method  employed  in  Interception  to  converge  to  the 
minimum  solution  of  Equation  (18).  The  technique  described  in  Section  4.1.1  for  address¬ 
ing  this  issue  has  been  employed  in  TurningCircle  .m.  It  was  found  by  trial-and-error  that 
setting  di  =  70000  and  (p  =  0.3  produced  faultless  results  for  0  ^  e  ^  0.9999,  regardless 
of  the  other  inputs.®  With  these  parameter  values,  TurningCircle  .m  will  terminate  after 
a  maximum  of  20  attempts  (10  starting  from  the  left  endpoint  plus  10  starting  from  the 
right  endpoint).® 

^See  http://www.wolfrcim.com/  for  more  information. 

^TurningCircle  .m  was  not  subjected  to  rigorous  testing.  However,  after  numerous  comparisons  with 
results  obtained  graphically  from  Equation  (18),  TurningCircle  .m  returned  no  false  outcomes. 

®To  obtain  an  indication  of  execution  time,  TurningCircle  .m  was  run  100  times  for  each  of  the  param¬ 
eter  values  used  to  generate  Figure  4(a)-(f)  (with  (5i  =  70000  and  tf)  =  0.3).  The  95%  confidence  interval 
for  the  mean  execution  time  in  CPU  seconds  is  [0.09,  0.40].  TurningCircle  .m  was  run  using  Mathematica 
6.0.2. 1  on  an  Apple  iMac  running  Mac  OS  X  10.4.11  with  a  2  GHz  PowerPC  G5  CPU  and  1  GB  of  RAM. 
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4.2.1  Examples  of  the  turning-radius  effect 

Six  examples  of  the  output  generated  by  TurningCircle  .m  are  displayed  in  Figure  4(a)- 
(f).  The  target’s  speed  increases  (relative  to  the  pursuer’s  speed)  from  Figure  4(a)  to 
Figure  4(f).  These  figures  illustrate  the  effect  of  turning  radius  as  they  encompass  all  of 
the  cases  discussed  in  Section  3. 

The  case  of  a  stationary  target  is  shown  in  Figure  4(a),  where  the  pursuer  can  turn 
either  left  or  right  to  intercept  the  target.  This  is  always  true  for  a  stationary  target, 
regardless  of  the  initial  heading  and  speed  of  the  pursuer,  provided  the  target  is  not  inside 
one  of  the  pursuer’s  turning-circles  (see  Proposition  8  of  Appendix  A. 2).  In  Figure  4(b) 
the  target  is  moving  slowly  enough  for  the  pursuer  to  intercept  the  target  before  it  enters 
either  turning-circle.  In  Figure  4(c)  Warning  1  (see  Section  4.1.3)  has  been  returned, 
and  consequently  the  pursuer  must  turn  away  from  the  target  to  perform  an  intercept. 
Warning  2  (see  Section  4.1.3)  has  been  returned  in  Figure  4(d)  and,  although  the  target 
has  greater  speed,  the  pursuer  must  still  turn  away  from  the  target  to  perform  an  intercept. 
In  Figures  4(e)  and  4(f),  the  target  has  sufficient  speed  for  the  pursuer  to  turn  towards  the 
target  to  perform  an  intercept  (the  solid  path),  however,  this  results  in  the  pursuer  being 
required  to  chase  the  target.  It  is  interesting  to  observe  that  increasing  the  relative  speed 
of  the  target  from  e  =  0.95  to  e  =  0.9999  leads  to  a  reduction  in  the  time  to  interception 
from  8.529  to  7.073,  respectively,  despite  the  pursuer  being  required  to  chase  the  target. 
Although  counterintuitive,  this  phenomena  occurs  because  as  the  target’s  speed  increases, 
the  time  taken  for  it  to  exit  the  turning-circle  decreases,  thus  enabling  an  intercept  to 
occur  earlier. 
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Solid  intercept  time:  2.5708 
Dashed  intercept  time:  8.35589 
Solid  intercept  point:  {2.,  2.} 
Dashed  intercept  point:  {2.,  2.} 
Solid  exit  point:  {2.,  1.} 

Dashed  exit  point:  {0.2,  -0.4} 


(a)  e  =  0,  which  corresponds  to  a  stationary 
target. 


Solid  intercept  time:  2.37461 
Dashed  intercept  time:  8.0044 
Solid  intercept  point:  {1.88127,  1.79435} 
Dashed  intercept  point:  {1.82281,  1.69311} 
Solid  exit  point:  {1.98636,  1.1646} 

Dashed  exit  point:  {0.209735,  -0.387235} 


(b)  e  =  0.1.  The  feasible  region  is  given  by  Case 
4  for  both  turns. 


Figure  Output  from  TurningCircl  e .  m,  which  is  an  implementation  of  the  Interception 
algorithm  as  a  Mathematica  package.  The  pursuer’s  path  to  interception  is  displayed  as  a 
solid  line  for  a  left  turn,  and  as  a  dashed  line  for  a  right  turn.  The  inputs  are:  Xin  =  (1,0), 
xljn  =  (2,  2),  Up  =  (1,  0),  and  ut  =  — (-^0.5/2,  ^1.5/2),  for  different  values  of  e  =  Ct/Cp 
(e  =  t)  and  e  =  O.lJ. 
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TurningCircle::warnl  :  Unable  to  perform  turn 
Dashed  intercept  time:  7.47982 
Dashed  intercept  point:  {1.55746,  1.23349} 
Dashed  exit  point:  (0.231286,  -0.360408} 


TurningCircle::warn2  :  Unable  to  perform  turn 
Dashed  intercept  time:  6.85649 
Dashed  intercept  point:  (1.23913,  0.682141} 
Dashed  exit  point:  (0.282434,  -0.30351} 


(c)  e  =  0.3.  The  feasible  region  is  given  by  Case  (d)  e  =  0.7.  The  feasible  region  is  given  by  Case 
3  for  the  left  turn,  and  Case  4  for  the  right  turn.  2  for  the  left  turn,  and  Case  4  for  the  right  turn. 

Figure  Output  from  TurningCircl  e .  m,  which  is  an  implementation  of  the  Interception 
algorithm  as  a  Mathematica  package.  The  pursuer’s  path  to  interception  is  displayed  as  a 
solid  line  for  a  left  turn,  and  as  a  dashed  line  for  a  right  turn.  The  inputs  are:  Xin  =  (1,0), 
=  (2,  2),  Up  =  (1,  0),  and  ut  =  — (  y^O.5/2,  y^1.5/2),  for  different  values  ofe  =  Cx/tp 
(e  =  0.3  and  e  =  0.7 ). 
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Solid  intercept  time:  8.52895 
Dashed  intercept  time:  6.63164 
Solid  intercept  point:  {-0.739756,  -2.7454} 
Dashed  intercept  point:  (1.12245,  0.480041} 
Solid  exit  point:  (0.0180475,  1.18913} 
Dashed  exit  point:  (0.31872,  -0.267977} 


(e)  e  =  0.95.  The  feasible  region  is  given  by  Case 
1  for  the  left  turn,  and  Case  3  for  the  right  turn. 


Solid  intercept  time:  7.07264 
Dashed  intercept  time:  6.59656 
Solid  intercept  point:  {0.0604049,  -1.35948} 
Dashed  intercept  point:  (1.10408,  0.44822} 
Solid  exit  point:  {0.000331144,  0.974267} 
Dashed  exit  point:  {0.326213,  -0.261074} 


(f)  e  =  0.9999.  The  feasible  region  is  given  by 
Case  1  for  the  left  turn,  and  Case  2  for  the  right 
turn. 


Figure  Output  from  TurningCircl  e .  m,  which  is  an  implementation  of  the  Interception 
algorithm  as  a  Mathematica  package.  The  pursuer’s  path  to  interception  is  displayed  as  a 
solid  line  for  a  left  turn,  and  as  a  dashed  line  for  a  right  turn.  The  inputs  are:  Xin  =  (1,0), 
=  (2,  2),  Up  =  (1,  0),  and  ut  =  — (■\/0.5/2,  y^l.5/2),  for  different  values  of  e  =  Ct/Cp 
(e  =  0.95  and  e  =  0.9999}. 
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5  Incorporating  classification  range 

In  maritime  surveillance  operations,  aircrew  can  classify  a  ship  if  the  aircraft  is  within  a 
distance  of  fd  from  the  ship,  where  fd  is  the  aircraft’s  classification  range  that  depends  on 
the  aircraft’s  sensor  suite,  weather,  sea  state,  et  cetera}^  In  this  section,  three  methods  are 
proposed  for  incorporating  classification  range  into  the  present  model:  an  exact  method, 
a  heuristic  method,  and  a  method  that  includes  an  angle  of  approach. 


5.1  Exact  method 

The  pursuer’s  classification  range  can  be  incorporated  into  the  present  model  by  transform¬ 
ing  the  interception  point  into  the  classification  point.  Let  Td  be  the  time  to  classification 
and 

XT(Td)  =  -h  eTdUT,  (30) 

the  position  of  the  target  at  the  time  of  classification.  The  classification  point  Xd  is  defined 
to  be  at  a  distance  of  fd  from  XT^{Td)  in  the  direction  of  Xout(Tc«)  —  xt(Tcj);  that  is,  the 
classification  point  results  from  the  most  direct  path  to  classification.  It  follows  that  Xd 
is  given  by 

Xd{Td)  =  xrr{Td)  +  |x  ^{T  i)- xt{T  i)\  ~  ^TC^d))  •  (31) 

To  determine  Td  it  is  necessary  to  replace  the  interception  point  x/  with  Xd  in  Equa¬ 
tions  (13)  and  (17).  As  a  consequence.  Equation  (17)  becomes  an  implicit  expression  for 
Xout(Td),  as  ti{Td)  and  Xd{Td)  are  also  functions  of  Xout(Td)-  Hence  the  following  system 
of  equations  must  be  solved  to  determine  the  minimum  time  to  classification: 

ti{Td)  =  \l\xd{Td)  -  Xcp  -  1,  (32) 

Xout(Td)  =  ^  tf{Td)xc  -h  {-l)'^ti{Td)  {xd{Td)  -  Xc)  x  k)  ,  (33) 

cos{Td  -  ti{Td))  =  (Xin  -  Xc)  •  (xout(Td)  -  Xc),  (34) 

where  Xd  is  given  by  Equation  (31). 

Since  Equation  (33)  is  an  implicit  expression  for  Xout(Tcz),  nearly  all  of  the  analysis 
in  Section  3  does  not  apply  to  Equations  (31)  to  (34),  and  so  the  Interception  algorithm 
cannot  be  easily  modified  to  accommodate  this  system.  Eurthermore,  the  analysis  of 
Equations  (31)  to  (34)  is  far  from  straightforward,  making  Constraints  (C1)-(C3)  difficult 
to  apply. 


5.2  Heuristic  method 

A  heuristic  method  that  accounts  for  the  pursuer’s  classification  range  and  returns  feasible 
times  to  classification,  can  be  obtained  by  initially  determining  the  minimum  time  to 

^°The  classification  range  in  this  section  refers  to  the  dimensionless  quantity  f^i  =  rcijrc,  where  is 
the  dimensional  classification  range. 
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interception  T  using  Equation  (29)  and  Interception.  An  approximate  classification  point 
is  then  derived  by  starting  at  x/(T)  and  moving  a  distance  of  fci  in  the  direction  of 
Xout(E)  —  x/(T).  The  resulting  classification  point  is  given  by 


Xci  =  < 


x/(r)  +  (xout(r)  -  x/(T))  , 


ti{T) 


rd  <  ti{T) 


rd  ^  ti{T), 


and  the  time  to  classification  is  simply^  ^ 


T-rd,  rd<ti{T) 
tc{T),  rd  ^  ti{T). 


(35) 


Observe  that  if  rd  ^  ti{T)  =  |xout(7")  ~x/(T)|  then  the  target  is  inside  the  pursuer’s 
classification  range  when  the  pursuer  exits  its  turn,  that  is,  at  time  tc(T). 


Equation  (35)  yields  a  feasible  time  to  classification.  However,  it  will  not  necessarily 
return  the  minimum  time  to  classification,  which  can  be  demonstrated  as  follows.  Since 
the  classification  point  is  closer  to  the  pursuer  than  the  interception  point,  the  pursuer’s 
classification  range  acts  as  if  to  increase  the  speed  of  the  pursuer.  Consequently  the 
pursuer  may  be  able  to  classify  the  target,  but  unable  to  perform  an  intercept.  In  this 
instance  the  heuristic  will  fail  because  the  approximate  classification  point  is  derived  from 
the  interception  point. 


5.3  Angle  of  approach 

The  pursuer’s  sensor  performance  or  its  path  can  usually  be  optimised  by  classifying  the 
target  at  a  specified  angle  of  approach  rj,  which  is  relative  to  the  target  [Mercer  et  al.  2008]; 
see  Eigure  5.  In  this  case,  the  classification  point  becomes^^ 

y^d{Td)  =  y^T{Td)  +  Id  (cos(r7),sin(r/),  0)  ,  (36) 

where  xt{Tci)  is  given  by  Equation  (30). 

To  determine  the  time  to  classification  Td,  it  is  necessary  to  replace  the  interception 
point  X/  with  x,,;  in  Equations  (13)  and  (17),  which  once  again  results  in  Equations  (32) 
to  (34).  Since  Equation  (36)  is  independent  of  Xout(7c0)  unlike  earlier,  the  analysis  in  Sec¬ 
tion  3  still  remains  valid  after  this  substitution.  This  is  because  evaluating  Equation  (32) 
using  Equation  (36)  yields  an  expression  for  ti{Td)  that  is  equivalent  to  Equation  (14), 
that  is, 

ii(Td)  =  \J {eTciY  +  -|-  d, 

where  all  that  has  changed  are  the  definitions  of  a  and  j3:^^ 

d  =  |x[},  -h  fci  (cos(?7),sin(r/),0)  -  Xc|^  -  1, 

^  =  ut  •  (x^  -h  fd  (cos(77),sin(?7),  0)  -  Xc)  . 

dimensional  units,  if  fcijC-p  <  ti{T)  then  Td  =  T  —  fci/C'p. 

^^This  representation  of  is  mathematically  consistent  with  the  technique  employed  in  Section  2.2. 

^^The  physical  meaning  of  a  and  j3  differs  from  that  of  a  and  (3. 


23 


DSTO-TR-2227 


Figure  5:  Two  feasible  paths  to  elassification  resulting  from  left  and  right  turns.  In  this 
figure,  the  pursuer  (P)  elassifies  the  target  (T)  at  7r/2  radians  on  the  target’s  left,  that  is, 
V  =  ^UT  +  in  Equation  (36). 


Hence  the  Interception  algorithm  can  be  easily  modified  to  accommodate  an  angle  of 
approach. 

The  angle  of  approach  can  be  chosen  relative  to  the  heading  of  the  target,  enabling 
classification  to  occur  on  the  left  or  right  of  the  target;  see  Figure  5.  Let  0ut  be  the 
heading  of  the  target  in  polar  coordinates  centred  on  the  target’s  position.  The  definition 
of  the  scalar  product  then  gives 


9 


ux 


arccos(uT  •  i),  ut  •  j  >  0 

27r  —  arccos(uT  •  i),  ut  •  j  ^  0, 


since  |ut|  =  1,  where  i  =  (1,  0,  0)  and  j  =  (0, 1,  0).  Now  set  rj  =  0ut  +  ^  in  Equation  (36). 
If  0  >  0  then  classification  will  occur  at  0  radians  to  the  left  of  the  target,  whereas  if  0  <  0 
then  classification  will  occur  at  9  radians  to  the  right  of  the  target. 
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6  Conclusion 

A  pursuer  intercepting  a  target  has  been  considered,  where  the  pursuer  and  target  move  at 
constant  speeds  in  two  dimensions  and  the  target’s  velocity  is  constant;  refer  to  Figure  1. 
The  pursuer  was  subjected  to  Constraints  (C1)-(C3);  in  particular,  the  pursuer  was  limited 
to  at  most  one  complete  turn.  The  minimum  feasible  path  to  interception  for  a  given 
turning  radius  was  sought. 

A  rigorous  analysis  of  Constraints  (C1)-(C3)  produced  the  Interception  algorithm  (re¬ 
fer  to  Figure  3)  that  can  be  used  to  systematically  search  the  feasible  region  for  the 
minimum  time  to  interception.  At  the  core  of  Interception  is  a  single  implicit  equation  for 
the  minimum  time  to  interception.  This  equation  is  valid  in  an  arbitrary  Cartesian  coor¬ 
dinate  system  and  encompasses  both  left  and  right  turns.  For  the  purpose  of  validation. 
Interception  has  been  implemented  as  a  Mathematica  package  that  displays  the  minimum 
feasible  path  to  interception,  as  shown  in  Figure  4(a)-(f). 

Three  methods  have  been  proposed  for  incorporating  classification  range  into  the 
present  model:  an  exact  method,  a  heuristic  method,  and  a  method  that  includes  an 
angle  of  approach.  These  methods  provide  simple  models  of  the  pursuer’s  sensor  perfor¬ 
mance.  The  Interception  algorithm  can  be  easily  modified  to  encompass  the  heuristic  and 
angle  of  approach  methods. 

The  point-to-point  and  unconstrained  movement  of  entities  in  some  simulations  of 
military  operations  is  unrealistic.  The  fidelity  of  this  representation  may  be  insufficient 
for  operations  research  studies.  The  Interception  algorithm  enables  the  effect  of  turning 
radius  to  be  incorporated  as  a  constraint  into  these  simulations,  improving  their  fidelity. 
Interception  is  also  straightforward  to  implement  when  compared  with,  for  example,  a  tra¬ 
ditional  flight  dynamics  model,  and  has  a  broad  range  of  applications  in  path  optimisation 
problems,  the  development  of  computer  games  and  robotics. 

To  summarise,  the  main  contributions  of  this  paper  are 

•  the  reduction  of  the  interception  scenario  to  a  single  equation; 

•  the  determination  of  the  feasible  region  as  a  function  of  the  scenario’s  inputs,  and; 

•  the  Interception  algorithm,  which  is  a  representation  of  an  entity’s  dynamics  in  terms 
of  its  turning  radius. 
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Appendix  A  Determining  feasible  times  to 
interception:  the  proofs 

The  primary  aim  of  this  appendix  is  to  rigorously  determine  which  times  satisfy  all  of 
Constraints  (C1)-(C3). 


A.l  Preliminary  results 

In  Section  2.3  it  is  claimed  that  the  solutions  to 


T-ti{T)  =  0,  (Al) 

T-ti{T)  =  TT,  (A2) 

are  given  by  Equations  (22)  to  (24).  These  solutions  were  constructed  by  solving 

T^  =  tf{T),  (A3) 

(r-7r)2  =  t2(T).  (A4) 


Solutions  of  Equations  (Al)  and  (A2)  are  also  solutions  of  Equations  (A3)  and  (A4), 
respectively,  however  the  reverse  implication  is  not  necessarily  true.  The  main  purpose  of 
the  following  Propositions  is  to  establish  when  solutions  of  Equations  (A3)  and  (A4)  are 
also  solutions  of  Equations  (Al)  and  (A2),  respectively. 


Proposition  1  Real  and  nonnegative  solutions  of  Equation  (A3)  also  solve  Equation  (Al). 


Proof  Let  T  be  a  real  and  nonnegative  solution  of  Equation  (A3).  This  implies  that  tf{T) 
is  also  real  and  nonnegative.  Therefore 

T=\T\  =  ^/^^=  =  \ti{T)\  =  ti{T). 


□ 


Proposition  2  Let  T  be  a  real  solution  of  Equation  (Af)  sueh  that  T  ^  tt.  Then 
T-TT  =  tl{T). 


Proof  Let  T  be  a  real  solution  of  Equation  (A4)  such  that  T  ^  tt.  This  implies  that  t‘f{T) 
is  also  real  and  nonnegative.  Therefore 

r  -  TT  =  IT  -  ttI  =  V(T  -  7r)2  =  =  \ti{T)\  =  tl{T). 


□ 
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A. 2  Feasible  times  to  interception  when  o  >  0 

Recall  that  the  target  is  initially  outside  the  pursuer’s  turning-circle  if  and  only  if  a  >  0. 
In  Section  2.3  it  was  shown  that  Constraint  (Cl)  is  satisfied  by  construction,  and  there 
are  essentially  two  cases  effecting  the  feasibility  of  times  to  interception  with  respect  to 
Constraints  (C2)  and  (C3)  when  a  >  0:  the  target  will  never  enter  the  interior  of  the 
pursuer’s  turning-circle  (e  =  0  or  /3  ^  0  or  ^  a),  and;  the  target  will  enter  the  interior 
of  the  pursuer’s  turning-circle  at  some  time  (0  <  e  <  1  and  /3  <  0  and  >  a). 

Proposition  3  Let  0  ^  e  <  1  and  a  >  0.  Then  Tq  as  defined  by  Equation  (22)  is  a  real, 
strietly  positive  number  and  uniquely  satisfies  tfiTfi)  =  0. 

Proof  The  conditions  of  the  proposition  directly  give  that  Tq  is  a  real  number  as  (1— 6^)0  > 
0.  If  /3  ^  0  then  Tq  >  0.  If  /?  <  0  then  the  inequality  —  \X\  +  -|-  Y  >  0  for  T  >  0 

yields  Tq  >  0.  It  can  be  shown  by  direct  substitution  that  Tq  solves  Equation  (A3)  and 
hence  tc(Eo)  =  0,  by  Proposition  1.  A  second  solution  to  Equation  (A3)  does  exist,  namely 


However,  by  using  the  aforementioned  method,  it  can  be  shown  that  Tq  ^  0.  □ 

Proposition  4  Let  a>0,  0<e<l,  /?<0  and  (3'^  >  a.  Then  Tq  ^  Tl- 

Proof  Recall  that  Tq  is  a  real  and  strictly  positive  number,  by  Proposition  3.  Under  the 
conditions  of  the  proposition, 

(1  ~  (eTo  T  fi)  =  e\J e^(/I^  —  cfi)  -\-  01  <,  -\-  \(3\  =0, 

and  so  Tq  <  —file.  It  follows  that  Tq  ^  Tl,  by  Proposition  3.  □ 

Proposition  5  Let  a  >  0  and  if  e  =  0  or  fi  0  or  fi'^  ^  a,  then  tfiT)  =  T  —  tfiT)  is  a 
strietly  inereasing  funetion  for  T  Tq.  If  0  <  e  <  1  and  fi  <0  and  >  a,  then  tfiT)  is 
strietly  inereasing  for  Tq  ^  T  ^  T^. 


Proof  Let  e  =  0  or  /3  ^  0  or  ^  a.  If  e  =  0  then  tfiT)  =  T—^/a  and  is  strictly  increasing. 
Now  let  0  <  e  <  1.  If  =  a  then  it  can  be  shown  that  t'fiT)  =  1  —  sgn(eT  -|-  fi)e  and 
hence  tfiT)  is  strictly  increasing.  Note  that 


,  ti{T)-e{eT  +  fi) 

’  ti{T) 

(A6) 

If  >  a  then  t'fiT)  =  0  iff 

ti{T)  =  e{eT  +  fi), 

(A7) 

which  has  the  solution 

T  — 

-(l-e2)/3+ V(l-e2)(/32-a) 

(A8) 

±  — 

e(l  -  e2) 
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If  /3^  <  a  then  Equation  (A8)  implies  that  tc{T)  has  no  real  stationary  points.  Let  ^  0 
and  >  a.  In  this  instance  it  can  be  shown  that  (1  —  e^)iT  —  Tq)  <  Tr  <  0  and  therefore 
tc{T)  has  no  stationary  points  for  T  ^  Tq.  Proposition  3  yields  t'^{To)  >  0,  which  concludes 
the  first  part  of  the  proof.  Now  let  0  <  e  <  1  and  /3  <  0  and  >  a.  Since  Tq  <  Tr  (by 
Proposition  4)  and  Tr  <  T,  tc{T)  has  no  stationary  points  for  Tq  ^  T  ^  Tr.  Then,  since 
t'PTo)  >  0,  tc{T)  is  strictly  increasing  for  Tq  ^  Tl.  □ 

Proposition  6  Let  a  >  0,  0  <  e  <  1,  /3  <  0,  >  a  and  Tr  ^  tt.  Then  feasible  times  to 

intereeption  T  oeeur  in  the  interval  Tq  ^T  ^Ti. 


Proof  Under  the  present  conditions,  ti{T)  is  real  and  non-negative  for  Tq  ^  T  ^  Tr  [see 
Equations  (14)  and  (21)]  and  therefore  Constraint  (C2)  is  satisfied.  By  Proposition  5,  tciT) 
is  strictly  increasing  for  Tq  ^  T  ^  Tr.  Purthermore,  tciTQ)  =  0  (by  Proposition  3)  and 
tc{TL)  =  Tl  ^  vr  (by  definition)  and  hence  Constraint  (C3)  is  satisfied  for  Tq  ^  T  ^  Tr.D 

Proposition  7  Let  a  >  0  and  e  =  0  or  /3  ^  0  or  ^  a.  Then  is  a  real  number 
that  uniquely  satisfies  tciTfi)  =  vr  and  feasible  times  to  intereeption  T  oeeur  in  the  interval 
Tq^t^  T^. 

Proof  If  e  =  0  then  T^-  =  TT  +  y/a  and  is  a  real  number.  If  /3  ^  0  then  tf{T)  is  an  increasing 
function  for  T  0  and  therefore 

tf{Tr)  +  -  a)  (efi)^  +  (l  -  e^)  a  >  0, 

implying  T-^  is  real.  If  ^  a  then  T-^  is  real  since 

^^(vr)  -|-  e^{0^  —  a)  =  (eir  P)  -|-  (l  —  {ct  —  /3^])  ^  0. 

Under  the  present  conditions,  tfiT)  is  real  and  non-negative  for  all  T  ^  0  [see  Equa¬ 
tions  (14)  and  (21)]  and  therefore  Constraint  (C2)  is  satisfied.  It  remains  to  show  that 
Constraint  (C3)  is  satisfied.  Under  the  conditions  of  the  proposition,  tfiT)  is  strictly  in¬ 
creasing  for  T  p  Tq  (by  Proposition  5)  and  since  tc(To)  =  0  (by  Proposition  3),  there  exists 
a  unique  solution  T  >  Tq  to  tfiT)  =  tt.  Recall  that  T  also  solves  Equation  (A4),  which 
has  two  real  solutions,  Tt^  and  T^,  that  are  given  by  Equations  (23)  and  (24),  respectively. 
It  follows  that  T  =  Ttt  or  T  =  T^.  Using  Equation  (A6)  it  can  be  shown  that  ffiTfi)  ^  0, 
however  this  violates  Proposition  5  and  hence  T  T^.  Therefore  tfiTfi)  =  tfiT)  =  vr. 
Pinally,  since  tciT)  is  strictly  increasing  for  T  ^  Tq,  Tq  <  T^r  and  0  ^  tciT)  ^  vr  for 
To  ^  T  ^  T^.  □ 


Proposition  8  Let  a  >  0  and  e  =  0  or  P  D  or  P'^  ^  a.  Then  T^^  is  a  real  number  and 
there  exists  a  feasible  solution  T  to  Equation  (18)  with  Tq  ^  T  ^  T^. 


Proof  Define  the  function  T(T)  to  be 

F{T)  =  cos(T  -  tl{T))  -  (Xin  -  Xc)  •  (xout(T)  -  Xc). 
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Note  that  a  zero  of  F(T)  is  a  solution  of  Equation  (18).  Under  the  conditions  of  the 
proposition,  ti{T)  is  real  and  continuous,  implying  F{T)  is  real  and  continuous.  Observe 
that  F{To)  ^  0  and  F{T.^)  ^  0,  since  by  definition 

-1  ^  (Xin  -  Xc)  •  (Xout(E)  -  Xc)  ^  1. 

It  follows  from  the  Intermediate  Value  Theorem  and  Propositions  3  and  7,  that  there  exists 
a  feasible  solution  T  to  Equation  (18)  with  Tq  ^  T  ^  T^.  □ 

Proposition  9  Let  a  >  0,  0  <  e  <  1,  P  <  0,  /3‘^  >  a  and  Tl  p  tt.  Then  T-n-  is  a  real 
number  that  uniquely  satisfies  tciTfi)  =  tt,  Tn-  and  there  exists  a  feasible  solution  T 

to  Equation  (18)  with  Tq  ^  T  ^  T^r. 

Proof  Rearranging  Tl  tt  leads  to  —{eir+P)  —  a,  implying  evr+Z?  <  0.  Therefore 

(evr  +  /3)^  —  a  >  {1  —  t^){P‘^  —  a),  that  is,  (ctt  +  /3)^  +  (1  —  e^)(a  —  /3^)  >  0  and  so 

Tjr  is  real.  Observe  that  tfiTo)  =  0  (by  Proposition  3),  tfiTL)  =  Tl  tt,  and  tc  is 
continuous  and  strictly  increasing  for  Tq  ^  T  ^  Tl  (by  Proposition  5).  It  follows  that 
there  exists  a  unique  solution  to  tfiT)  =  tt  with  Tq  <  T  ^  Tl-  Recall  that  T  also  solves 
Equation  (A4),  which  has  two  real  solutions,  Tn-  and  that  are  given  by  Equations  (23) 
and  (24),  respectively.  It  follows  that  T  =  or  T  =  T^.  Using  Equation  (A6)  it  can  be 
shown  that  t'JfT-,^)  ^  0,  however  this  violates  Proposition  5  and  hence  T  T^.  Therefore 
tciTn)  =  tc{T)  =  TT,  Tl  Ttt  and  0  ^  tfiT)  ^  tt  for  Tq  ^  T  ^  T^^.  Then  by  following 
the  reasoning  in  the  proofs  of  Propositions  7  and  8,  it  can  be  shown  that  there  exists  a 
feasible  solution  T  to  Equation  (18)  with  Tq  ^  T  ^  T^-.  □ 

Proposition  10  Let  a  >  0,  0  <  e  <  1,  /3  <  0,  >  a  and  Tji  ^  tt.  Then  is  a  real 

number  that  satisfies  tc(T7r)  =  tt  and  feasible  times  to  intereeption  T  oeeur  in  Tq  ^  T  ^  Tl 
and  Tfj  ^  T  ^  T^r. 


Proof  Rearranging  Tr  ^  tt  leads  to  eTT  +  P  \/ P‘^  —  a,  implying  evr  +  /3  >  0.  As  in  the 
proof  of  Proposition  9,  it  follows  that  T^r  and  T^^  are  real.  Since  Tl  <  Tr  ^  tt,  feasible 
times  to  interception  occur  in  Tq  ^  T  ^  Tr,  by  Proposition  6.  Recall  that  T  is  the  unique 
stationary  point  of  tfiT)  and  T  >  Tji  [see  Equations  (21)  and  (A8)].  Eurthermore 


tc{T) 


e^(/3^  —  a) 

tfiT) 


and  hence  tfiT)  has  a  local  minimum  at  T.  Observe  that 


tc{T) 


-/3+  \/(l  -  e^){P‘^  -  qQ 

e 


(A9) 


(AlO) 


and  0  <  tc{T)  <Tr  ^  tt,  implying  that  Constraint  (C3)  is  satisfied  for  Tr  ^  T  ^  T.  Since 
0  <  tc{T)  <  TT  and  tc  is^continuous  and  strictly  increasing  for  T  >  T,  there  exists  a  unique 
solution  T  >  T  to  tc(T)  =  tt.  By  following  the  argument  in  the  proof  of  Proposition  9, 
it  can  be  shown  that  T  =  T^r.  Eurthermore,  Tn-  >  T  and  Constraint  (C3)  is  satisfied  for 

T  <  r  ^  T^.  □ 
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Proposition  11  Let  be  a  eomplex  number  or  en  +  (3  ^  0.  Furthermore,  let  a  >  0, 
0<e<l,  /3<0  and  >  a.  Then  <  vr  <  Tr  and  feasible  times  to  intereeption  T 
only  oeeur  in  Tq  ^  T  ^  T^. 

Proof  By  Propositions  9  and  10,  if  ^  vr  or  Tr  ^  vr  then  T,r  is  real,  which  under 
the  present  conditions  implies  that  Tr  <  vr  <  Tr.  It  follows  from  Proposition  6  that 
feasible  times  to  interception  occur  in  Tq  ^  T  ^  Tr.  It  remains  to  prove  there  does 
not  exist  any  feasible  times  to  interception  for  T  >  Tr.  Constraint  (C2)  is  not  satisfied 
for  times  Tl  <  T  <  Tr.  Now  let  T  ^  Tr.  If  T,r  is  complex  then  Equation  (23)  yields 
|e7r  +  P\  <  Y^(l  —  e^){(3‘^  —  a).  It  follows  that  — (evr  +  (3)  +  y^(1  —  e^){(3‘^  —  a)  >  0,  as 
X  ^  |a;|  for  any  real  number  x.  Then  from  Equation  (AlO), 

e{tc{T)  -  vr)  =  -(evr  +  /?)  +  v^(l  -  e‘^){/3‘^  -  a)  >  0,  (All) 

and  so  tc(T)  >  vr.  Recall  that  tc(T)  is  the  global  minimum  of  tc{T)  for  T  >  Tr  and 
tc{TR)  =  Tr  >  TT,  therefore  Constraint  (C3)  is  violated  for  T  ^  Tr.  If  evr  +  /3  ^  0  then 
Equation  (All)  implies  that  Constraint  (C3)  is  violated  for  T  ^  Tr.  □ 

Proposition  12  Let  T^r  be  a  real  number,  en  +  P  >  0,  a  >  0,  0  <  e  <  1,  (3  <  0,  (3"^  >  a 
and  Tr  <  TT  <  Tr.  Then  tciT-P)  =  tc{TT^)  =  vr  and  feasible  times  to  intereeption  T  oeeur 
in  To  ^  T  ^  Tr  and  T^  ^  T  ^  T^^. 

Proof  Since  Tr  <  n,  feasible  times  to  interception  occur  in  To  ^  T  ^  Tr,  by  Proposition  6. 
Since  T^r  is  real.  Equation  (23)  yields  |e7r  +  /3|  ^  y^(1  —  —  a).  As  evr  +  /3  >  0,  it 

follows  that 

-(evr  +  /3)  +  \/{l  -  e^)(/3^  -  a)  ^  0, 

and  so  tc{T)  ^  vr  [see  Equation  (All)].  Recall  that  tc{T)  is  the  global  minimum  of  tc{T) 
for  T  >  Tr.  Eurthermore,  tciTR)  =  Tr  >  vr  and  Tr  <  T.  Hence  there  exist  real  solutions 
Ti,  T2  to  tc{Ti)  =  TT  (for  i  =  1,2)  such  that  Tr  <  Ti  ^  T  and  T  ^  T2.  Undei^the  present 
conditions.  Equation  (A4)  has  exactly  two  real  solutions  T^-,  T^r,  and  Ti,  T2  also  solve 
Equation  (A4).  Therefore  T^-  =  Ti  and  T^^  =  T2,  and  feasible  times  to  interception  occur 
in  T^  ^  T  ^  T,r-  n 

A. 3  Feasible  times  to  interception  when  o  ^  0 

Recall  that  the  target  is  initially  inside  (or  on  the  boundary  of)  the  pursuer’s  turning-circle 
if  and  only  if  a  ^  0.  Then  Constraint  (C2)  is  only  satisfied  for  times  T  ^  Tr,  because  the 
target  does  not  leave  the  pursuer’s  turning-circle  until  T  =  Tr. 

Proposition  13  Let  a  ^  0  and  0  <  e  <  1.  Then  Tr  ^  0  and  Tr  p  0. 

Proof  The  proof  follows  from  the  definitions  of  Tr  and  Tr  and  the  inequality  |/3|  ^ 
V/32  -  a.  □ 

Proposition  14  Let  a  ^  0  and  0  <  e  <  1.  Furthermore,  let  T^r  be  a  eomplex  number  or 
evr  -|-  /3  ^  0.  Then  feasible  times  to  intereeption  do  not  exist. 
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Proof  Since  a  ^  0  and  0  <  e  <  1,  feasible  times  to  interception  cannot  exist  in  the 
domain  T  <  Tr  (see  Proposition  13).  By  following  the  proof  of  Proposition  11,  it  can  be 
seen  that  feasible  times  to  interception  do  not  exist  for  T  ^  T^.  □ 

Proposition  15  Let  a  ^  0,  0  <  e  <  1,  Tq  be  a  real  number  and  (3  0.  Then  Tq  ^  0  and 

satisfies  tciTfi)  =  0,  and  tc{T)  ^  0. 

Proof  The  conditions  of  the  proposition  and  the  definition  of  Tq  give  Tq  0  and  tfiT^)  =  0 
(by  Proposition  1).  Since  Tq  is  a  real  number,  —  a)  +  afi  t).  This  implies 

-  |/3|  +  V(1  -  e2)(/32  -  a)  ^  0. 

Recall  that  tfiT)  has  a  unique  minimum  at  T  >  Tr;  see  Equation  (A8).  The  above  inequal¬ 
ity  in  conjunction  with  the  present  conditions  and  Equation  (AlO)  yield  min'r>R^  A(E)  = 

tc{T)  ^0.  □ 

Proposition  16  Let  a  ^  0  and  0  <  e  <  1.  Furthermore,  let  Tq  be  a  eomplex  number  or 
P  <  0.  Then  miuTyTn  tc{T)  >  0. 


Proof  Let  /3  <  0.  Then  Equation  (AlO)  gives  min'r>'r^  A(T)  =  tc{T)  >  0.  Now  let  Tq  be 
a  complex  number.  Then  following  the  proof  of  Proposition  15  leads  to 

-fi  +  \/(l  -  e2)(/32  _  a)  -\P\  +  \/(l  -  e2)(/32  _  a)  >  0, 

and  therefore  minR^R^  tc(E)  =  tfiT)  >  0,  by  Equation  (AlO).  □ 

Proposition  17  Let  a  ^  0,  0  <  e  <  1,  Tr  ^  vr,  Tq  and  be  real  numbers  and 

P  0.  Then  tfiTo)  =  tfiTo)  =  0,  tciTfi)  =  tt,  feasible  times  to  intereeption  T  oeeur 

in  Tfi  ^  T  ^  Tq,  and  there  exists  a  feasible  solution  to  Equation  (18)  with  Tq  ^T  ^  T,r- 

Proof  Recall  that  for  T  >  Tr,  tc  is  real,  continuous  and  has  a  unique  minimum  of  tfiT). 
In  addition,  tfiTfi)  =  Tr  ^  0  and  tfiT)  ^  0,  by  Propositions  13  and  15.  Hence  there  exist 
real  solutions  Ti,  T2  to  tfiTp  =  0  (for  i  =  1,  2)  such  that  Tr  ^  Ti  ^  T  and  T  ^  T2.  Under 
the  present  conditions.  Equation  (A3)  has  exactly  two  real  nonnegative  solutions  Tq  and  Tq 
[see  Equation  (A5)],  and  it  follows  that  Tq  =  Ti  and  Tq  =  T2,  by  Proposition  1.  Likewise, 
it  can  be  shown  that  tciTfi)  =  tt  and  Tq  <  T^r.  It  can  be  shown  using  Equation  (A6) 
that  tc  is  strictly  decreasing  for  Tr  ^  T  ^  Tq  and  strictly  increasing  for  Tq  ^  T  ^  T^^. 
Then,  since  0  ^  tfiTfi)  ^  tt,  feasible  times  to  interception  occur  in  Tr  ^  T  ^  Tq  and 
To  ^  T  ^  T,r-  Observe  that  tfiT)  <  0  for  Tq  <  T  <  Tq  and  therefore  Constraint  (C3)  is 
violated  in  this  interval.  The  existence  of  a  solution  to  Equation  (18)  in  Tq  ^  T  ^  T^-  can 
be  seen  by  following  the  proof  of  Proposition  8.  □ 

Proposition  18  Let  a  ^  0,  0  <  e  <  1,  Tr  >  tt,  Tq  and  T^^  be  real  numbers  and  P  0. 

Then  tc(To)  =  tc(To)  =  0,  tciTfi)  =  tcfiTfi)  =  tt,  feasible  times  to  intereeption  T  oeeur  in 

Ttt  ^  T  ^  Tq,  and  there  exists  a  feasible  solution  to  Equation  (18)  with  Tq  ^  T  ^  T^-. 
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Proof  The  proof  is  the  same  as  that  of  Proposition  17,  with  the  following  addition.  Since 
tciTji)  >  TT  and  tc(7"o)  =  0  (refer  to  the  proof  of  Proposition  17),  it  follows  that  tc(7",r)  = 
with  Tr  <  <  Tq.  Furthermore,  Constraint  (C3)  is  violated  for  Tr  ^  T  <  and 

feasible  times  to  interception  occur  in  ^  T  ^  Tq.  □ 

Proposition  19  Let  a  ^  0,  0  <  e  <  1,  Tr  ^  tt,  T-,^  be  a  real  number  and  eir  +  (3  >  0. 

Furthermore,  let  Tq  be  a  eomplex  number  or  (5  <t).  Then  tciT^^)  =  vr  and  feasible  times  to 

intereeption  T  oeeur  in  Tji  ^  T^r. 

Proof  Recall  that  feasible  times  to  interception  do  not  exist  for  T  <  Tr,  and  observe  that 
0  <  tc{T)  <  Tr  ^  TT  (by  Proposition  16).  The  proof  can  be  completed  by  following  the 
proof  of  Proposition  10.  □ 

Proposition  20  Let  a  ^  0,  0  <  e  <  1,  Tr  >  tt,  T-,^  be  a  real  number  and  eir  +  (3  >  0. 

Furthermore,  let  Tq  be  a  eomplex  number  or  (3  <t).  Then  fc(F,r)  =  tc{T-K)  =  and  feasible 

times  to  intereeption  T  oeeur  in  T-,^  ^T  ^  T^-. 

Proof  Recall  that  feasible  times  to  interception  do  not  exist  for  T  <  Tr,  and  observe 
that  tc{T)  >  0  (by  Proposition  16).  The  proof  can  be  completed  by  following  the  proof  of 
Proposition  12.  □ 
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Appendix  B  An  explanatory  flow  chart  of  the 

Interception  algorithm 

Below  are  descriptions  of  the  errors,  warnings  and  parameters  that  appear  in  Figure  Bl(a)- 
(b): 

Error  1  is  returned  if  the  target’s  speed  is  greater  than  or  equal  to  the  speed  of  the 
pursuer. 

Error  2  is  returned  if  the  root-finding  method  used  to  solve  Equation  (18)  fails  to  con¬ 
verge,  when  a  solution  is  known  to  exist;  see  Propositions  8  and  9  from  Appendix  A. 2. 

Warning  1  is  returned  if  the  root-finding  method  used  to  solve  Equation  (18)  fails  to 
converge.  In  this  instance,  the  pursuer  is  unable  to  intercept  the  target  before  it 
enters  the  turning-circle  (if  a  >  0),  and  the  target  is  moving  too  slowly  to  leave  the 
turning-circle  before  the  pursuer  has  performed  a  complete  turn;  see  Case  3  from 
Section  3. 

Warning  2  is  returned  if  the  root-finding  method  used  to  solve  Equation  (18)  fails  to 
converge.  In  this  instance,  the  pursuer  is  unable  to  intercept  the  target  before  it 
enters  the  turning-circle  (if  a  >  0),  and  the  target  is  moving  too  slowly  to  leave  the 
turning-circle  in  time  for  the  pursuer  to  attempt  an  intercept;  see  Cases  1  and  2 
from  Section  3. 

To  is  the  time  to  interception  if  the  pursuer  does  not  turn  [see  Equation  (22)]. 

Tl  is  the  time  for  the  target  to  enter  the  pursuer’s  turning-circle  [see  Equation  (21)]. 

Tr  is  the  time  for  the  target  to  exit  the  pursuer’s  turning-circle  [see  Equation  (21)]. 

Tjr  is  the  time  to  interception  if  the  pursuer  takes  tt  units  of  time  to  turn  [see  Equa¬ 
tion  (23)]. 

Ttt  is  the  time  to  interception  if  the  pursuer  takes  vr  units  of  time  to  turn  [see  Equa¬ 
tion  (24)]. 
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Input:  direction  of  the  turn, 
speeds  of  the  pursuer  &  target, 
initial  positions  of  the  pursuer  &  target, 
initial  velocity  of  the  pursuer,  and 
the  velocity  of  the  target 


No- 


Continued 

on 

Figure  Cl  (a) 


(  ^  ) 


(a)  The  cases  where  the  pursuer  is  able  to  perform  a  complete  turn  without  the  target  entering 
the  pursuer’s  turning-circle.  The  other  cases  are  presented  in  Figure  Bl(b). 


Figure  Bl:  The  Interception  algorithm  (continued  on  Figure  Bl(b))  for  determining  fea¬ 
sible  times  to  interception  T,  described  using  words  rather  than  symbols;  refer  to  Figure  3 
36  for  the  symbolic  version.  Descriptions  of  the  errors,  warnings  and  parameters  can  be  found 
on  page  35. 
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(b)  The  cases  where  the  target  will  enter  the  pursuer’s  turning-circle  before  the  pursuer  completes 
a  turn.  The  other  cases  are  presented  in  Figure  Bl(a). 


Figure  Bl:  The  Interception  algorithm  (continued  from  Figure  Bl(a))  for  determining 
feasible  times  to  interception  T,  described  using  words  rather  than  symbols;  refer  to  Fig¬ 
ure  3  for  the  symbolic  version.  Descriptions  of  the  errors,  warnings  and  parameters  can 
be  found  on  page  35. 
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Appendix  C  The  Interception  algorithm  when 

o  ^  0 

The  Interception  algorithm  for  the  case  when  the  target  is  initially  inside  (or  on  the 
boundary  of)  the  pursuer’s  turning-circle  (a  ^  0)  is  presented  in  Figure  Cl(a)-(b).  Refer 
to  Appendix  A. 3  for  the  derivation  of  Interception  for  the  a  ^  0  case. 

Below  is  a  brief  description  of  the  algorithm  for  this  case: 

(A) :  Is  the  target  stationary  and  initially  on  the  boundary  of  the  pursuer’s  turning-circle? 

(B) :  Will  the  target  remain  inside  the  pursuer’s  turning-circle  after  the  pursuer  has 
performed  a  complete  turn? 

(C) :  Is  the  target  near  the  boundary  of  the  pursuer’s  turning-circle  and  heading  out  of 
the  turning-circle? 

(D) :  Can  the  pursuer  feasibly  attempt  an  intercept  immediately  after  the  target  leaves 
the  pursuer’s  turning-circle? 

(E) :  In  this  case,  the  target  is  near  the  boundary  of  the  pursuer’s  turning-circle  and 
heading  out  of  the  turning-circle.  If  the  pursuer  has  not  intercepted  the  target 
immediately  after  the  target  leaves  the  turning-circle,  then  the  target  will  be  far 
enough  away  from  the  turning-circle  for  the  pursuer  to  feasibly  perform  a  complete 
turn  before  attempting  an  intercept. 
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(a)  The  cases  where  the  target  is  stationary  and  on  the  boundary  of  the  pursuer’s  turning-circle,  or 
where  the  target  is  not  near  the  boundary  of  the  turning-circle  or  heading  into  the  turning-circle. 
The  other  cases  are  presented  in  Figure  Cl(b). 


Figure  Cl:  The  continuation  of  the  Interception  algorithm  for  determining  feasible  times 
to  interception  T,  where  the  initial  position  of  the  target  is  inside  (or  on  the  boundary  of) 
the  turning- circle  of  the  pursuer.  The  Errors  and  Warnings  are  described  in  Section  4.1.3, 
and  refer  to  the  Notation  section  on  page  xi  for  descriptions  of  the  parameters. 
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Yes 


I 

(  ^  ) 

(b)  The  case  where  the  target  is  near  the  boundary  of  the  pursuer’s  turning-circle  and  heading  out 
of  the  turning-circle.  The  other  cases  are  presented  in  Figure  Cl(a). 

Figure  Cl:  The  continuation  of  the  Interception  algorithm  for  determining  feasible  times 
to  interception  T,  where  the  initial  position  of  the  target  is  inside  (or  on  the  boundary  of) 
the  turning- circle  of  the  pursuer.  The  Errors  and  Warnings  are  described  in  Section  4-T3, 
and  refer  to  the  Notation  section  on  page  xi  for  descriptions  of  the  parameters. 
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Appendix  D  The  Mathematica  package: 

TurningCircle . m 


BeginPackage ["TurningCircle ‘ "] 

TurningCircle : :usage=" TurningCircle [eps ,purinit ,tarinit ,purvel,tarvel] 

The  dimensionless  parameter  eps=ct/cp: 

-  cp:  pursuer's  speed 

-  ct:  target’s  speed. 

It  is  assumed  that  0  <=  eps  <  1 . 

TurningCircle  also  requires  the  following  inputs  to  be  defined  as 
2-vectors : 

-  purinit :  pursuer’s  position  at  t=0 

-  tarinit :  target’s  position  at  t=0 

-  purvel:  pursuer’s  velocity  at  t=0  ( | purvel | ==1) 

-  tarvel:  target’s  velocity  ( | tarvel | ==1) . 

All  variables  have  been  scaled: 

X  by  r,  t  by  r/cp,  purvel  by  cp,  and  tarvel  by  ct;  here  r=radius 
of  the  pursuer’s  turning-circle. 

An  optional  argument  BlackAndDashed->True  will  display  the  paths 
as  black  and  dashed  curves . " 


TurningCircle :: errl="eps  must  satisfy  0  <=  eps  <  1" 
TurningCircle :: err2=" I purvel I  and  |tarvel|  must  be  1" 
TurningCircle: :err4=" I purinit-tarinit I  must  be  >  0" 
TurningCircle :: err5= "Unable  to  compute  intercept  times" 
TurningCircle :: err7="FindRoot  failed  to  converge" 
TurningCircle : :warnl= "Unable  to  perform  turn" 
TurningCircle : :warn2= "Unable  to  perform  turn" 

Options  [TurningCircle] ={BlackAndDashed->False} 

Begin ["TurningCircle ‘Private ‘ "] 
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i={l,0,0};  j={0,l,0};  k={0,0,l>; 

(^Centre  of  the  two  turning-circles,  n=0,l*) 
xc [n_] : =xin+(-l) * (n+1) *Cross [up,k] 

(^Pursuer’s  intercept  points,  where  T=tC+tL  is  the  total  intercept  time, 
eps=cs/ ca*) 

xI[T_,eps_] : =xs+eps*T*us 

(*Time  from  turning-circle  exit  point  to  intercept  point,  tL*) 

A [n_] : =If [Abs [Norm [xs-xc [n] ] -1] >NumTol ,Dot [xs-xc [n] ,xs-xc  [n] ] -1 ,0] 

B  [n_] : =Dot [us , xs-xc [n] ] 

tL  [T_, eps_ ,n_]  :=Sqrt  [(eps*T)  ''2+2*eps*B  [n]  *T+A  [n]  ] 

(^Pursuer’s  turning-circle  exit  points*) 
xout [T_ , eps_ , n_] : = 

(xl [T, eps]  +tL [T, eps ,n] *2*xc [n] +(-l) "n^tL [T,eps ,n] *Cross  [xl  [T, eps] -xc  [n] ,k] )/ 
(l+tL[T,eps,n] "2) 

(*Value  of  T  corresponding  to  tC=0*) 

TZero  [eps_  ,n_]  :  =  (eps*B  [n]  +Sqrt  [(eps*B  [n] )  ''2+(l-eps''2)  *A  [n]  ] )  /  (l-eps*2) 
TZeroBar  [eps_,n_]  :  =  (eps*B  [n]  -Sqrt  [(eps*B  [n] )  *2+(l-eps~2) *A [n] ] ) / (l-eps''2) 

(*Value  of  T  corresponding  to  tC=Pi*) 

TPi [m_ , eps_ , n_] : = 

(eps*B  [n]  +Pi+  (-1)  ~  (m+1)  *Sqrt  [  (eps*Pi+B  [n] )  ''2+  (l-eps*2)  *  (A  [n]  -B  [n]  ''2)] )  / 
(l-eps''2) 

(*The  left  and  right  endpoints  of  the  solution’s  domain  when  B[n]  <  0  AND 
B[n]''2  >=  A[n]  AND  0  <  eps  <  1*) 

TLeftRight [m_, eps_ ,n_] :=(-B [n] +(-l) " (m+1) *Sqrt [B [n] *2-A [n] ] ) /eps 

(^Pursuer’s  interception  times:  T=tC+tL*) 

Interception [eps_ ,n_] :=( 

(*Check  Assumptions*) 

Which [ 

eps<0 1  I eps>=l , 

Message [TurningCircle : : errl] ;  Abort  [] , 

Norm [up] !=1| | Norm [us] !=1, 

Message [TurningCircle : : err2] ;  Abort  [] , 

Norm [xin-xs] <NumTol , 

Message [TurningCircle : : err4] ;  Abort  [] 

]; 

If  [A[n]>0, 

AlphaPosInterception [eps ,n]  , 
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AlphaNeglnterception [eps ,n] 

] 


(*If  Target  starts  outside  the  turning-circle*) 

AlphaPosInterception [eps_ ,n_] :=Module [ 

(*Local  Variables*) 

{T , TLef t , TRight} , 

If  [0<eps<l&&B  [n] <0&&B [n] ~2>A [n] , 

(*The  harder  case*) 

If [TLef tRight  [0 , eps ,n] <Pi , 

TRight=TLef tRight [0 , eps , n]  , 

TRight=TPi [1 , eps ,n] 

]; 

(*Look  for  solution  TZero  <=  T  <=  T-  OR  TZero  <=  T  <=  TPi+*) 
T=SearchForT [TZero [eps ,n] , TRight , eps ,n] ; 

If  [NumberQ [T] , 

Return [T] , 

If  [TLef tRight [0 , eps , n] >=Pi , 

Return [Message [TurningCircle : : err7]  ] 

]; 

]; 

(*If  the  Target  will  not  leave  the  turning-circle  in  time  ...  *) 

If  [  (eps*Pi+B  [n] )  *2+  ( l-eps''2)  *  (A  [n]  -B  [n]  ''2)  <0  I  I  eps*Pi+B  [n]  <=0 , 

Return [Message [TurningCircle : : warnl] ] 

]; 

(*If  a  solution  has  not  been  found  AND  an  Abort  []  has  not  occurred, 
look  for  solution  T+  <=  T  <=  TPi+  DR  TPi-  <=  T  <=  TPi+*) 

If [TLef tRight [1 , eps ,n] >Pi , 

TLef t=TPi [0 , eps , n] , 

TLef t=TLef tRight [1 , eps ,n] 

]; 

T=SearchForT [TLef t , TPi [1 , eps ,n] , eps ,n] ; 

If  [NumberQ [T] , 

Return [T] , 

Message [TurningCircle : : warn2] 


43 


DSTO-TR-2227 


]; 


(*The  easy  case*) 

(*Look  for  solution  TZero  <=  T  <=  TPi+*) 
T=SearchForT [TZero [eps ,n] ,TPi [1 , eps ,n] , eps ,n] ; 

If  [NumberQ [T] , 

Return [T] , 

Return [Message [TurningCircle : : err7]  ] 

]; 


] 

(*End  AlphaPosInterception*) 

] 

(*If  Target  starts  inside  the  turning-circle*) 

AlphaNegInterception [eps_ ,n_] :=Module [ 

(*Local  Variables*) 

{T,TLeft}, 

(*If  Target  is  stationary  and  on  the  rim  of  the  turning-circle*) 

If  [eps==0&&A  [n] ==0 , 

Return  [ArcCos [Dot [xin-xc [n] , xs-xc [n]  ]  ]  ] 

]; 

(*If  Target  will  not  leave  turning-circle*) 

If [ (eps*Pi+B [n] ) *2+ ( l-eps~2) * (A [n] -B  [n] ~2) <0 1  | eps*Pi+B  [n] <=0 , 

Return  [Message [TurningCircle : : warnl]  ] 

]; 

(*If  Target  is  not  near  the  rim  or  heading  into  the  turning-circle*) 
If  [  (eps*B  [n] )  •'2+  (l-eps*2)  *A  [n]  <0  |  |  B  [n]  <0 , 

If [TLef tRight [1 , eps ,n] >Pi , 

TLef t=TPi [0 , eps , n] , 

TLef t=TLef tRight [1 , eps ,n] 

]; 

T=SearchForT [TLeft ,TPi [1 , eps ,n] , eps ,n] ; 

If  [NumberQ [T] , 

Return [T] , 

Return [Message [TurningCircle : : warn2]  ] 
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]; 


(*Else  if  Target  is  near  the  rim  and  heading  out  of  the  turning-circle*) 
If [TLef tRight [1 , eps ,n] >Pi , 

TLef t=TPi [0 , eps , n] , 

TLef t=TLef tRight [1 , eps ,n] 

]; 

T=SearchForT [TLef t jTZeroBar [eps ,n] ,eps,n] ; 

If  [NumberQ [T] , 

Return [T] 

]; 

T=SearchForT [TZero [eps ,n] ,TPi [1 , eps ,n] , eps ,n] ; 

If  [NumberQ [T] , 

Return [T] , 

Return [Message [TurningCircle : : err7]  ] 

]; 


] 

(*End  AlphaNegInterception*) 

] 

(^Finding  the  solution*) 

SearchForT [xL_ ,xR_ , eps_ ,n_] : =Module [ 

(*Local  Variables*) 

{Fraclnit=70000 , FracPer=0 . 3 , Frac , FR, T} , 

If [xR<NumTol , Return [] ] ; 

Frac=FracInit ; 

While  [ ! VectorQ  [FR] &&Frac>l , 

FR=CheckAbort  [ 

FindRoot [Cos [T-tL [T,eps ,n] ] ==Dot [xin-xc [n] ,xout [T,eps ,n] -xc  [n] ] , 
{T, (xR+(Frac-l) *xL) /Frac} , 

EvaluationMonitor : >If [Im[T] !=0| |T<xL| | T>xR, Abort []] ] , 

Null] ; 

Frac=FracPer*Frac ; 

]; 

(*If  unsuccessful,  try  again  from  the  right*) 

If  [! VectorQ [FR] , 
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Frac=FracInit ; 

While [ ! VectorQ  [FR] &&Frac>l , 

FR=CheckAbort [ 

F indRoot  [Cos  [T-tL  [T , eps , n] ] ==Dot  [xin-xc  [n] , xout  [T , eps , n] -xc  [n] ] , 

{T, (xL+(Frac-l)*xR)/Frac}, 

EvaluationMonitor :  >If  [Ini[T]  !  =0  |  |  T<xL  |  |  T>xR,  Abort  []  ]  ]  , 

Null] ; 

Frac=FracPer*Frac ; 

]; 

]; 

If  [VectorQ [FR] ,T/.FR] 

(*End  SearchForT*) 

] 

(*Time  from  turning-circle  entry  point  to  turning  circle  exit  point,  tC*) 
tC[T_,eps_,n_] :=If  [ 

Dot [xout [T , eps , n] -xin , up] >0 , 

T-tL  [T, eps, n] , 

2*Pi- (T-tL [T , eps , n] ) 

] 

(^Pursuer’s  intercept  times:  T=tC+tL,  correcting  for  the  "2*Pi  effect"*) 
TotalTime [T_,eps_ ,n_] :=N [tC [T,eps ,n] +tL [T,eps ,n] ] 

(*The  arcs  of  the  pursuer’s  turning  circles*) 

PathArc [T_ , eps_ , n_] : =Module [{Thin} , 

If [Dot  [xin-xc  [n] , j ] >0 , 

ThIn=ArcCos [Dot [xin-xc [n] , i] ] , 

ThIn=2*Pi-ArcCos [Dot [xin-xc [n] , i] ] 

]; 

If [EvenQ [n] , 

{Thin , ThIn+tC  [T , eps , n] } , 

{ThIn-tC[T,eps,n] ,ThIn} 

] 

] 

(*Define  graphics  objects*) 

Circles [T_ , eps_ ,n_ , colour_] := 

Graphics [{Thick, colour , Circle  [{xc [n]  [[1]] ,xc[n]  [ [2] ] } , 1 ,PathArc [T,eps ,n] ] }] 
Pursuer Initial : = 

Graphics [{Thick, Arrowheads [Medium] , 

Arrow  [{{xin  [  [1]  ]  -up  [  [1]  ]  ,  xin  [  [2]  ]  -up  [  [2]  ]  } ,  {xin  [  [1]  ]  ,  xin  [  [2]  ] »]  }] 
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InterceptPoint [T_ , eps_ ,n_ , colour_] : = 

Graphics [{Thick , colour , Arrowheads [Medium] , 

Arrow [{{xout[T,eps,n] [[!]] ,xout [T,eps,n] [[2]]},{xl  [T,eps]  [[!]] ,xI[T,eps]  [[2]]}}]}] 

TargetHeading [T_ , eps_] := 

Graphics [{Thick, Arrowheads [Medium] , 

Arrow [{{xs [[!]] ,xs [[2]]},{xl [T,eps] [[!]] ,xI[T,eps]  [[2]]}}]}] 

(*The  possible  path(s)  and  time(s)  to  interception:  eps=cs/ca*) 

TurningCircle [eps_ , acinit_ , shinit_ , acvel_ , shvel_ , OptionsPattern [] ] : =Module  [ 

(*Local  Variables*) 

{Disp , DutputStyle , PathStyle , TO , Tl} , 

(*0ptions*) 

If [OptionValue  [BlackAndDashed] ==True , 

OutputStyle={"Solid" , "Dashed"} ; 

PathStyle={Black, Dashed}; , 

(*Def ault*) 

OutputStyle={"Blue" , "Red"} ; 

PathStyle={Blue , Red} ; 

]; 

(*Pursuer’s  entry  point  at  t=0.  Global  Variable.*) 
xin=Append [acinit , 0]  ; 

(*Target’s  position  at  t=0.  Global  Variable.*) 
xs=Append [shinit , 0] ; 

(*Pursuer’s  velocity  at  t=0.  Note:  |up|==l.  Global  Variable.*) 
up=Append [acvel , 0] ; 

(*Target’s  velocity.  Note:  |us|==l.  Global  Variable.*) 
us=Append [shvel , 0] ; 

(*Numerical  Tolerance.  Global  Variable.*) 

NumTol=10~(-6) ; 

T0=lnterception[eps,0] ; 

Tl=Interception  [eps , 1] ; 

(*Generate  graphics*) 

Which [ 

! NumberQ [TO] &&NumberQ [Tl]  , 

Print  [DutputStyle [ [2] ], "  intercept  time:  " ,TotalTime  [Tl , eps , 1] ] ; 

Print  [DutputStyle [ [2] ], "  intercept  point:  " ,N  [{xl  [Tl ,eps]  [  [1] ] ,xl  [Tl , eps]  [  [2] ] }] ] ; 
Print [DutputStyle [ [2] ], "  exit  point:  " ,N [{xout [Tl , eps , 1] [ [1] ] ,xout  [Tl , eps , 1]  [ [2] ] }] ] 
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Disp=Show [InterceptPoint [T1 , eps , 1 ,PathStyle [ [2] ] ] ,TargetHeading [T1 , eps] , 

Circles [T1 , eps , 1 ,PathStyle [ [2] ] ] ,PursuerInitial , 

Graphics [Circle [{xc [0] [  [1] ] ,xc  [0]  [  [2] ]}]]];, 

NumberQ  [TO] && ! NumberQ  [Tl] , 

Print [DutputStyle [ [1] ], "  intercept  time:  " ,TotalTime  [TO, eps , 0] ] ; 

Print  [DutputStyle [ [1] ], "  intercept  point:  " ,N  [{xl  [TO , eps] [ [1] ] ,xl [TO , eps] [ [2] ] }] ] ; 
Print [DutputStyle [ [1] ], "  exit  point:  " ,N[{xout [T0,eps,0] [[1]] ,xout  [T0,eps,0]  [[2]]}]] 
Disp=Show [InterceptPoint [TO, eps , 0,PathStyle [ [1] ] ] ,TargetHeading [TO, eps] , 

Circles  [TO , eps ,0 ,PathStyle [ [1] ] ] ,PursuerInitial , 

Graphics [Circle [{xc [1] [ [1] ] ,xc [1] [ [2] ] }]  ]  ]  ; , 

NumberQ [TO] &&NumberQ [Tl] , 

Print  [DutputStyle [ [1] ], "  intercept  time:  " ,TotalTime  [TO, eps , 0] ] ; 

Print [DutputStyle [ [2] ], "  intercept  time:  " ,TotalTime  [Tl , eps , 1] ] ; 

Print  [DutputStyle [ [1] ], "  intercept  point:  " ,N  [{xl  [TO , eps]  [  [1] ] ,xl  [TO , eps]  [  [2] ] }] ] ; 
Print [DutputStyle [ [2] ], "  intercept  point:  " ,N [{xl  [Tl , eps]  [ [1] ] ,xl  [Tl , eps]  [ [2] ] }] ] ; 
Print [DutputStyle [ [1] ], "  exit  point:  " ,N[{xout [T0,eps,0] [[1]] ,xout [T0,eps,0] [[2]]}]] 
Print  [DutputStyle [ [2] ], "  exit  point:  " ,N  [{xout  [Tl ,eps , 1]  [  [1] ] ,xout  [Tl , eps , 1]  [  [2] ] }] ] 
Disp=Show [InterceptPoint [TO, eps , 0,PathStyle [ [1] ] ]  , 

InterceptPoint [Tl , eps , 1 ,PathStyle [ [2] ] ] ,TargetHeading [TO, eps] , 

TargetHeading [Tl , eps] , Circles [TO , eps , 0 , PathStyle [ [1] ] ] , 

Circles  [Tl , eps , 1 , PathStyle [ [2] ] ] , Pursuerinitial] ; , 

! NumberQ [TO] && ! NumberQ [Tl]  , 

Message [TurningCircle : : err5] ;  Abort  [] ; 

]; 

(^Display  graphics*) 

Disp 

] 

End  [] 

EndPackage  [] 
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Appendix  E  The  points  Xjn  and  Xout  in  polar 
coordinates  centred  on  Xc 


It  may  be  useful  to  express  the  points  Xin  and  Xout  in  polar  coordinates  centred  on  Xc, 
that  is,  find  values  of  Oia  and  0out  such  that 

Xin  =  Xc(n)  +  (cos(6»in(n)),sin(0in(n))), 

Xout(r,n)  =  Xc(n)  +  (cos(6'out(r,n)),sin(0out(r,n))), 

where  9^  and  0out  depend  on  which  direction  the  pursuer  turns  (n  =  0  or  n  =  1),  and  T 
is  the  minimum  time  to  interception  returned  by  Interception;  see  Figure  3.  Values  for 
0in  and  0out  can  be  calculated  using  the  arctan  function,  however  this  may  not  take  into 
account  which  quadrants  of  Xc(n)  the  points  Xin  and  Xout  are  in.  Alternatively,  expressions 
for  Bin  and  0out  can  be  deduced  from  Figure  1  such  that  the  resulting  points  are  in  the 
correct  quadrant.  These  expressions  are 


since  |xin 


Bm{n) 


-  Xc(n)| 


arccos((xin  -  Xc(n))  •  i),  (xjn  -  Xc(n))  •  j  >  0 

27r  -  arccos((xin  -  Xc(n))  •  i),  (xjn  -  Xc(n))  •  j  ^  0, 
1,  where  i  =  (1,  0, 0)  and  j  =  (0, 1,  0),  and 
9ont{T,n)  =  6-M  +  i-irUf), 


where  tc  is  given  by  Equation  (28). 

By  definition,  the  direction  of  increasing  polar  angle  is  anticlockwise.  Therefore,  since 
the  pursuer  moves  in  an  anticlockwise  direction  on  one  turning-circle  and  in  a  clockwise 
direction  on  the  other  turning-circle,  it  can  be  shown  that 


6'in(0)  ^  6»out(F,  0),  n  =  0 
9out{T,l)  ^  Bin{l),  n  =  l. 
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